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Abstract 

The purpose of the PQIC Optimal Surface Mitigated Multiple Targeting System (PQIC 
OSMMTS) is to define the equipment and processing necessary to produce, in real time, 
an error- bounded, self-monitoring and self-adjusting, likelihood- based Target Position 
Report for arbitrarily many self-identifying targets in a two-dimensional grid. Each 
target sends identifying information to an array of sensors strategically placed in its 
vicinity to maximize the likelihood that the system will calculate a position report as 
accurately and precisely as possible. 
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1. Mission Statement 

The purpose of the PQIC™ Optimal Surface Mitigated Multiple Targeting System (PQIC OS- 
MMTS 1 ) is to define the equipment and processing necessary to produce, in real time, an error-bounded, self- 
monitoring and self-adjusting, likelihood-based Target Position Report for arbitrarily many self-identifying 
targets in a two-dimensional grid 2 . Each target sends identifying information to an array of sensors strate- 
gically placed in its vicinity to maximize the likelihood that the system will calculate a position report as 
accurately and precisely as possible. The OSMMTS uses analytical and ad-hoc mitigation and optimization 
techniques to reduce the error bounds on the target report to a practical minimum. This memorandum 
provides the baseline technical documentation of the analytical methods, construct guidelines, quantifi- 
cation methods, mitigation and optimization techniques, and programming details for implementing the 
PQIC OSMMTS. Additional supporting methods, documented separately, that supplement, complement, 
and refine the analytical methods found herein shall be issued in an ongoing series through PQIC technical 
documentation. 



0 Process Quality Improvement Consulting, PQIC, Optimal Surface Mitigated Multiple Targeting System, and OSMMTS 
are trademarks of PQI Consulting, P. O. Box 425616, Cambridge, MA 02142-0012 USA, info@pqic.com. 

1 While the formal name of the system is PQIC OSMMTS for copyright, trademark, service mark, and patent purposes, 
the name most often used in this, and subsequent technical, marketing, and promotional documentation shall be the OSMMTS 
System, or simply the OSMMTS. 

2 A separate, yet related, targeting system for a three-dimensional grid will be covered in a separate set of PQIC documen- 
tation. 
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2. The OSMMTS System Interface 

The OSMMTS System Interface consists of one Principal Application Specific Integrated Circuit (PASIC) 
Central Processing Unit (CPU), generically referred to as the PCPU, a set (at least four) of remote sensing 
SDU's that send information to the PCPU, and a database of statically stored data that the PCPU accesses 
for parameter 3 data, algorithm exceptions, and other information, which are used to produce the Target 
Position Report (TPR) 4 , as well as supporting reports as the implementation determines (see Figure 2). The 
PCPU, SDU's, and any database systems must be coordinated on and agree with an absolutely maintained 
time system, accurate to at least twice the precision of the anticipated Target Position Report. In this 
respect, the accuracy of the Target Position Report depends more on the maintenance of the time system 
than it does on the algorithmic introduction of error due to calculation round-off or lack of precision. 

A Target Position Report is generated whenever a SDU sends a coordinated stream 5 of timing infor- 
mation to the PCPU. Since different SDU will send information at slightly different times about the same 
target, an absolute timing schedule must be used to ensure valid comparison of timing data from the SDU 
set. 

A Target T may only initiate a signal to the SDU set when t = 0 modf, where f = 10 n /p cycles in a 
10 n Hz PCPU, where there are p signals per second For example, for a 1 GHz PCPU, if a target sends a 
signal to the SDU set every \ second, then p = 2, and £ = ^ = 1Q tf 10 2 = 10 9 - lo Sio 2 . 

The Effective Range of the OSMMTS System is the maximum time for this receive/query /confirm 
period. It measures the farthest a target may be away from the closest qualifying set of SDU's and still be 
detected by the system. 

A complete Signal Period, i.e., £ = 10 n /p cycles in a 10 n Hz PCPU, consists of six Phases, each 
encompassing an interaction between the PCPU, the SDU set, and the parameter database (see Figure 1). 

1. Receive, during which the PCPU receives the detected signal information from the SDU set. This phase 
must last as long as the effective range, plus overhead time for communications between the SDU set 
and the PCPU. The information passed during this phase consists of: 

(a) SDU ID 

(b) Target ID 

(c) Time Of Signal Detection 

The SDU and Target ID are static codes used throughout all phases and signal periods. If either 
the SDU ID or the Target ID changes during a signal period, it must be through a formal change 
management process incorporated into the particular implementation of the OSMMTS System. It shall 
be the responsibility of the OSMMTS System implementation to ensure that changing SDU ID and/or 
Target ID are linked properly for inference purposes. 

The Time of Signal Detection is relative to the common absolute timing mechanisms in the OSMMTS 
System. 

2. Query, during which the PCPU queries the sending SDU for a confirmation code to ensure communica- 
tion integrity. If the confirmation code sent by the SDU is not correct (see the next phase), the PCPU 
queries the SDU again for the proper confirmation code. This is repeated up to a tunable number of 
iterations. If no correct confirmation code is received in the allotted time, the SDU is deactivated. 

3. Confirm, during which the PCPU receives and processes the confirmation code sent by the SDU. It is 
during this phase that any required re-transmissions are also requested, received, and disposed. 

4. Process, during which all calculations are completed to produce the Target Position Report, and sub- 
sequent reports for evaluation, quantification, and adjustment purposes. 

5. Report, during which the Target Position Report and supporting information are made available on 
output channels, and during which any auxiliary communications with the SDU's is completed. 

6. Sync, during which no processing activity is scheduled. This is useful when coordinated processing 
activities require synchronized signal periods. 

One signal period begins when the previous one ends. The sync phase may be used to coordinate 
any overhead processing issues to implement this requirement. 

3 Parameter data consists of Demerit, History, and Confirmation values for each SDU. See Section 45 on page 38 for more 
information. 

4 See Section 42 on page 37. 

5 In the OSMMTS, "coordinated" means "timed to a cycle coordinated by a master clock." 
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Signal Timing 




Figure 1: The OSMMTS Signal Timing Cycle 



A special SDU, the Calibrating SDU (CSDU), is placed where the PCPU can always receive its signal 
as clearly as possible, i.e., without corruption. This CSDU calibrates the system by being a known distance 
from all other SDU's, and allows for a consistent, guaranteed error-free TPR by which all other SDU's axe 
evaluated for accuracy and consistency. 
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Figure 2: The OSMMTS Processing System 



3. Target Position Calculations 

4. Introduction 

Suppose there were two Surface 6 Detection Units (SDU) labeled I and II that detect a signal that propagates 
through a uniformly dense environmental medium 7 encompassing the SDU and a Target T 8 . The signal 
from the target is detected at times t\ and £2, where the detection is made under a common methodology 9 . 
Suppose further that the SDU are a known distance d > 0 apart, and that T and the SDU are si > 0 and 
S2 > 0 distances apart, respectively. Finally, let 9 be the angle (measured positively) between T and the 
SDU (see Figure 5). In the case that T is on the line segment between I and II, take 9 - | , i.e., the triangle 
formed by I, II, and T is degenerate, with two 0 angles between I and T, and II and T. If I, II, and T are 
otherwise collinear, then take 0 = 0; the other angles will be irrelevant. 

The calculations that follow will depend on axioms that collectively assert the OSMMTS System. These 
axioms formally state the conditions or assumptions under which the OSMMTS System may be implemented. 

Axiom 1 Given d, and the positions of I and II, the only observed data, and therefore the only basis for 
any calculations found in the OSMMTS, are the times t± and t 2 of signal detection. 



6 The detection units are called "surface" units because all coordinates under the OSMMTS may be thought of as (x,y } z) € 
R 3 , were 2 = 0. Even after rotations and translations, the coordinates may be thought of as planes or surfaces embedded in 
R 3 . 

7 The term "environmental medium" is intentionally generic. The analytical and implementation aspects of the OSMMTS 
apply equally well to signals broadcast through the air, under water, in gas-filled containers, in vacuums, or in any other 
medium. 

However, the medium will be assumed to be of uniform density to ensure the validity of the formula: 

Distance = (Linear Speed) * Time 

The "uniformly dense" restriction may be relaxed in subsequent OSMMTS System versions. 

8 Even though this development applies to a single Target T, all analytical methods in the OSMMTS may be extended to 
any number of individual targets with distinguishable signals. 

9 The "common methodology" refers to the electrical, mechanical, physical, or other well-defined, reproducible, and deter- 
ministic means of declaring a detection complete at an SDU. Whatever particular methods are used, the implementation must 
be uniform and constantly applied to each SDU in the OSMMTS System to ensure the integrity of algorithmic calculation. 
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5. Derivation 

A critical element of the OSMMTS is to calculate a Target Position Report, which contains a position report 
based on the arrived times at the SDU's. Consider arrival times t\ and t 2 at two SDU, named I and II, for 
the same Target T. 



II 



• T 



Figure 3: Relative placement of I, II, and T 



For the triangle formed by I, II, and T depicted in Figure 5, the Law of Cosines gives 

d 2 = s\ + s\ - 2s x s 2 cos 8 (X) 

Let / = 1 1 — £2 j c = the speed with which the signal linearly propagates through the environmental 
medium, and let t be the time at which the Target T sends the signal, with units in common. Note how 
ti = t only when T is identically collocated with either I or II. The term "speed" refers to the instantaneous 
change in distance per unit time, and it is expressed in terms of two-dimensional Euclidean geometry. 

Axiom 2 For the purposes of the OSMMTS, it will be assumed, without loss of generality, that U > t, for 
i = 1,2. 



Then 



or 



Hence, 



fc = (*i - * 2 ) c 

= ((*!-*)- (*2 ~ t)) C 

= (*i - t) c - (f 2 - t) c 

5?-2s 1 S 2 + 5l = (/c) 2 

5 i + S 2 ~" 2siS2 cosfl - d 2 = 0 
s\ + si - 2sis 2 - (fc) 2 =0 
2s 1 s 2 (l -cos0) =cP -{fcf 



or 



S\S 2 = 



# - (fcf 
2(1 - cos 6) 



= K 



(2) 
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since 9 ^ 2n7r, for any integer n. If I, II, and T are collinear, e.g., the Target T is on the line between I 
and II, then take 9 = n (as in the case of a degenerate triangle). Assume otherwise 10 that I, II, and T are 
not collinear. 
Note how 

d*-(/c) 2 rf» - (/c) 2 _ Q m 
/C ~2(l-cos^) =* cos6 '- 1 Yk (3) 

Here, 

k > 0, when d > |/| c 
k = 0, when d = j/j c 
k < 0, when d < |/| c 
since d and c are both necessarily positive. 
Now s\ — $2 = fc gives us 

si (s x — fc) — K 
s i — f cs i — « = 0 

(-$)'- (*)' 



Hence, 



si = ^ ( /c ± \/4k + (/c) 2 ) and s 2 = 2<C 

V 7 fc±yj4K + (f C y 



(4) 



The choice of sign in (4) depends on the signs of / and k; whichever makes $i > 0 and s 2 > 0 should 
be used. 

However, since si > 0 and 52 > 0, then k ^ 0, hence 

2k = ( 2k \ ( fc*y/AK+(fcf 

fc±y/AK+{fc)' 2 \fc±y/4K+(fcf ) [fcTy/^Hfc) 2 

2k \ f fcTy/^K+jfcf 

fc±y/4 K+{fc)* J \fcT\/4K+Uc) 2 

2*(fcT\/4K+(fc) 2 ) 
~ -4k 

(/c Tx /4k+(/c) 2 ) 

_ 2 



So 



5l = i (/c± ^4k+(/c) 2 ) and s 2 = i (-/c ± ^4k + (/c) 2 ) 



(5) 



Expression (5) shows that the are functions of ac, /, and c. Since k is a function of d, /, c, and 0, 
we have the Si as functions of d, /, c, and 9. The values of d and c are fixed by the circumstances and by 
nature, and the / value is observed; this means the s» may be thought of as functions of 9 given d and the 
data t\ and £ 2 - 

We have 

4« + (/c) 2 = 4(^^) + (/c) 2 

= o ( d 2 -Uc? \ , (l-cosfl)(/c) 2 

M l-cos0 y 1-COS0 
_ 2d 2 -2(/c) 2 + (l-cos fl)(/c) 2 

1— cos0 
= 2d 2 -(l-fcosfl)(/c) 2 

~ 1— COS0 

hence 



10 If I, II, and T are collinear, yet the Target T is not on the line between I and II, then while it would make sense to take 
6 = 0, this would make k ill-defined. For the purposes of the OSMMTS development, assume I, II, and T are collinear only 
when the Target T is on the line between I and II. 
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Note how 



Hence, 



i (/c ± ^4/c + (/c) 2 ) -j(- /c ± + (/c) 2 ) 

(/c) 2 ± 2/c^4/ c+(/c) 2 + 4k + (/c) 2 
-(-/c) 2 ± 2/c/u+(/c) 2 -4k- (/c) 2 
±/ CN /4« + (/c) 2 



Sl + «2 = 



*1 

«l-«2 



_ ±/<V4*+(/c) : 



(6) 



»/>i»yi) 



5i \ ^JJ (x 2 ,ift) 

(z, y) 



c 2 



Figure 4: Intersection of Circles C\ and C2 



Since 5i and S2 are distances from T to I and II, respectively, then T must be at the intersection of 
two circles centered at I and II (see Figure 4). In coordinate geometry, suppose I is at (xi,yi), II is at 
(22,2/2), and T is at (x,y). The origin for the coordinate system is arbitrary, and its position may be chosen 
for convenience in any particular- implementation. Then 

(x-xi) 2 + (y-yi) 2 = $i 
(x - x 2 ) 2 + (y - 2/2) 2 = 5 2 



must be satisfied simultaneously for a given 0. Hence, from (6), 
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x 2 - 2xxi + x\ + y 2 - 2yyi + y\ = s\ 

x 2 - 2xx 2 ± x\ + y 2 - 2yy 2 + y| = a| 

2x (x 2 - xi) + 2y (y 2 - yi) = (a? - s|) + [x\ - x\) + (yf - y 2 ) 

2x (x 2 - xi ) + 2y (y 2 - y x ) = ±/c^4«+ (/c) 2 + (x?, - xf ) + (y| - y 2 ) 
Now if yx = y 2 , then we have 

2x (x 2 - xi) = ±fc\]±K+{fcf + {x\ - x\) 



or 



_ ±fc^4K + (fc) 2 + (x 2 -x 2 ) 
2{x 2 - xi) 

which is a vertical line. If yi ^ y 2 , then we have 

2x (x 2 - xi) + 2y (y 2 - y x ) = ±fcyfeit + (/c) 2 + (x 2 , - x 2 ) + (y| - y 2 ) 



or 



= - x ( X *- X A 

\V2-y\) 



+ 



±/ C ^4k+(/c) 2 + (xl - x?) + (y 2 - yj) 



2(y2 - yi) 



which is a line with slope - anc * y-intercept ±/c ^ 4K+(/c ^y 2 ^ l ^ 1 ^^ 2 To summarize, let 



Li : 



±fcsjAK+Uc? + (xl-x\) 

x = 2 (x 2 -x!) ; f 1 = ^ 2 

±fcy/*K+(fc) 2 + (xl-x*)+(yl-y*) 
2(y 2 -yi) 



(7) 



is a line in x and y. Hence, as 6 varies, the set of points of intersection between Circles I and II is a line 
defined by L\ (see Figure 5). 



A 0*1, yi) 




Figure 5: Foci of Points as 9 Varies 

If a third SDU, named III, is added to the circumstance, at position (x 3 ,y 3 ), with sides relabeled to 
account for the third SDU, and which is not collinear to the positions of SDU I and II (see Figure 6), and 
subject to the axioms mentioned earlier, e.g., see Axiom 2, then another line of points of intersection L 2 , 
defined as in (7), may be expressed in terms of the same equations as before, this time for Circles I and III 
(see Figure 7). In particular, the equation for Li becomes 



x = ■ ■ — * 



Li : 



2(X2-Xt 



y = - X ( + 



dr/iCv/4/ci+(/ic) 2 + (x^z;)-h(y|-y?) 
2(y 2 -yi) 



yi = 2/2 

2/1 7^2/2 



t 
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and 



x 



L 2 : 



±/2cy r 4/c 2 + (/ 2 c) a -f(xl-x;) 
~ 2(xs-xil ~_J 

±/2C % /4/C2 + (/2C) 2 + (x^X^) + (v|-V?) 

2(y 3 -yi) 



v=-*(s) + 

for distinct K\ and /C2, themselves functions of distinct 0i and #2- 



2/i = 2/3 

2/1 7*2/3 



/// 



S22" - - - yr> 



Figure 6: Third SDU Added 

However, since Su = S12 (see the common side in Figure 6), then by (5), we have a link between ki 
and K2, namely 



4/ci + (ficj 



or 



Therefore, we have 



= 1 f-/ 2 c ± ^4js2 + (/ 2 c) 2 ) 
(/1 + /a) c ± + (Ac) 2 = ±^4k 2 + (/ 2 c) : 



±/>c N /4« l +(/ic) 2 +(xl-i?) 

X = - 



r J 2(X2— xi 

1 1 *=-*tet) + 



2(y2-yi) 



2/1 = 2/2 

2/1 7*2/2 



(8a) 



and 



x 



_ / 2 c((/i-f/ 2 )c±V4 M +(/ic) 2 )-h(xl-x?) 

_ 2(13-11) > 

/2c((/i+/2)c±V4/c 1 4-(/ic) 2 )+(xi-x;)+(yg- y ;) 

2(y 3 -yi) 



2/1 = 2/3 
> 2/1 7* 2/3 



(9a) 



y = -x + 

y \ya-yi y 

Claim 3 // I, II, and III are not collinear, then L\ and L<i are not parallel. 

Proof. Let fiu be the slope of the line between SDU i and j, e.g., /Z23 = x 2 It 3 , or fi23 = +00 for a 
vertical line. Since I, II, and III are not collinear, we have 

M12 / ^23 / ^13 

with allowances for dealing with infinities, and £1,2:2, X3 are not all the same value. 

Suppose 2/1 7* 2/2 and yi ^ y 3 . Then the slope of L\ is - ( j^ffj* ) > 311(1 f° r ^2 it is - ( f^Eff ) • If 



_ / x 2 - Xj \ _ _( x 3 - Xi \ 

V 2/2 - yi / V 2/3 - 2/1 / 
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Figure 7: Intersection of Foci Locates T 



then either 

2/2 - 2/1 2/3 - 2/1 



X 2 — X\ £3 — X\ 

or 

2/3 — 2/i 

+00 = , when X\ — £ 2 

or 

2/2-2/1 , , 

= +00, when x\ = x 3 

X2 — X\ 

If V~V — v J~Z l > tnen M12 = Vi3\ a contradiction. 

If +00 = jj^Ijj^ , then xi = x 3 , since 1/1 7^ 2/3. This is also a contradiction, since X\ = x 2 under this 
situation. 

If l\Zx\ = + 00 , tnen £1 = #2) since yi 7^ 1/2, yet another contradiction, since x\ = x$ under this 
situation. 

If either y\ = y 2 and y 2 ^ 2/3 3 or y 2 = 2/3 and yi 7^ y 2 , then the slope of L\ would be +00, and for L 2 it 
would be - ( yglyl )> or ^ ne °ther wa Y around. In either case, these slopes could not be equal in any sense. 

Finally, if both yi = y 2 and y\ =2/3, then I, II, and III would be collinear, since they would all have 
the same y- value; this is also a contradiction. ■ 

By Claim 3, let 6q be the value of 0 for the point of intersection between L\ and L 2 . When yi ^ y 2 
and 2/1 ^ y 3 , let 



771 1 



f x 2 -xi \ fxz-xA 
— _ and m 2 = - 

\ 2/2 - 2/1 / V 2/3 - 2/1 / 



with 



= ±/iC^4* 0 + (Ac) 2 4- (j = xl) + (gj - y 2 ) 

1 2{y 2 - yi ) 

= ±f 2 cy/4Ko + (/ 2 c) 2 + (*j - j) + (^ - y 2 ) 

2 2(y 3 -y x ) 



for «o corresponding to #o 
Then 



L\ : y = m\x + 61 
L 2 :y = m 2 x + 62 
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Claim 4 When I, II, and III are not collinear, and y x ^ y 2 and y\ ^2/3, the point of intersection between 
Lx andL 2 is (x 0) y 0 ) = (;£=fe.»»i (^E&j) +6,). 
Proof. We have 

miXo + b\ = 7712^0 + &2 



or 



and 



fr 2 ~ h 

771 x — 7712 



J/oLj = mix 0 -h6i 

= (mi-m 2 + m 2 )(^5^)+<, 1 

= ^-6i+6i+m 2 (^) 

= soli, 

■ 

Claim 5 When I, II, and III are not collinear, and y\ = j/ 2 onrf i/i ^ J/3 ( tne point o/ intersection between 
L\ and L 2 is 



(a;o,yo) = 



V 

Proo/. We have 



2(y 3 -yi) 



V 



±flCy/4Ko+(f lC ) 2 + (zl-xl) 
2(x 2 -xi) 



\ 



he ((/i +f 2 )c± yjiio + (/ic) 2 ) + (x§ - x?) + (y| - y\) 



Li : x = 



and 



/x 3 - Xi \ 
L 2 :y = -x[ 

V V3 - Vi J 



+ 



±/icy4/co + (/ic) 2 + (xj-xf) 

2(x 2 -xi) 

/ 2 c ((/i + / 2 ) c + ^4/c 0 + (/ic) 2 ) + (xl - x\) + (2/| - »?) 

2(2/3 -yi) 



Hence 
2>o|/, a 



/2c((/ 1 +/2)c±V4i C o-K/ic) 2 )+(x^-x?)+(yj-y?) 

2(y 3 -yi) 



c%/4/cq + (/ic) 2 -{-(x^-x' 



2(ya-yi) 



xa- 




1 + 


"/2c((/i+/ 2 )c± N /4/c 0 +(/ic) 2 ) + (x 2 -x 2 )+(y 2 - 


-*?) 


ya - 


-yi > 


2(y 3 -yi) 




x a - 


-Xi> 


1 + 


■/2c((/i+/ 2 )c± x /4« 0 + (/ic) 2 ) + (x 2 -x 2 ) + (y 2 - 


-»?) 


x 2 - 




2(y 3 -yi) 



2(W3-»i) 



/ 2 c ((/i + h) c ± ^/4/co + (/ic) 2 ) + (x§ - x 2 ) + (y 3 2 - j/ 2 ) ^ 
- (i/ic^/co + (Ac) 2 + 



Claim 6 PV7ien J, XT, and J/J are not collinear, and yi ^ 2/2 and 2/1 =2/3, £ne point of intersection between 
L\ and L 2 is 

I /2C 2 (/i+/ 2 )± N /4k 0 +(/ic) 2 +(x 2 -x 2 ) ^ 

' 2(x 3 ^xT l ' 

x ±ficyJtK Q + { jrf + {xl- x\) + (2/| - y?) 

V { - (/2C ((/l + /,) C ± ^ + (/ lC ) 2 ) + (Xi - X?)) (g=£) y y 



(zo,2/o) = 



14 
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Proof. We have 



and 



Li :y = -x 



X 2 — X\ 

2/2 ~ 2/1 



+ 



±ficy/^p + (Ac) 2 + (j - x\) + (»J - »?) 



2(2/2-2/i) 



Hence 



/ 2 c ((/1 + / a ) c ± ^4k 0 + (fief) + (a§ - x 2 ) 

L 2 : x = ^ 

2(x 3 -xi) 



±/icV4ico+(/ lC ) 2 +(zg- a ;)+(v;-y;) 

2(y 2 -t/i) 



_ _ ^ /2c((/i+/ 2 )c±V4ko+(/ic) 2 ) + (x|-^) 



2(x 3 -xi) 



/ \V2-yi / 



/ /2c((/i+/ 2 )c±V4K 0 +(/ 1 c) 2 ) + (x 2 -> a; 2 ) \ X 

y 2(2/2-2/1) J \x3-x1 y 



+ 



+ 



jb/ 1 cV4/c 0 + (/ lC ) 2 +(x 2 -x 2 ) + (y 2 -y 2 ) 
2(2/2-2/1) 

±/icV4K 0 +(/ lC ) 2 +(x 2 -x 2 )+(yi-y 2 ) 
2(2/2-1/1) 



2(2/2 -yi) 



±hcy]±K 0 + ( / lC ) 2 + (x?> - x?) + (y| - 1/ 2 ) 

" (f2C ((/! + / 2 ) C ± ^4* + (/1C) 2 ) + (Xl - X 2 )) (gEjl) 



The fourth possibility consistent with Claims 4-6, namely yi = y 2 and y\ =2/3, has been excluded by 
the assumption that I, II, and III are not collinear. 

Note how in all Claims 4-6, (x 0 ,yo) is a function of constants and the Sij (via the which are in turn 
functions of 0. This means 0 O may be expressed implicitly through the point (xo,2/o)» 

If 

Vi = (x u y x ) - (x 0 ,2/o) 

and 

V 2 = (x 2 ,y 2 ) - (x 0 ,2/o) 

then I Vi | 2 = 5n and \V2\2 = $21, so that 

Vi - V 2 = |Vi| 2 I Va | 2 cos 0 O = 5iis 2 i cos 0 O 
This proves the following claim. 

Claim 7 0 O = Arc cos f 

\«11*21 / 

The quadrant of 0 O must be considered when calculating the inverse cosine in Claim 7. Once established, 
we then have 

S11 (0o) S21 (Oo) cos 0 O = (xi - x Q (0 o )) (x 2 - x 0 (0 O )) + (2/1 - 2/0 (0o)) (2/2 - 2/o (0o)) (10) 
However, from (2), 

su (0o) S21 (0o) = ko 

and, from (3), 

2«o 



so 



Sll (0o) 521 (0o) COS 0 O = « 0 ( 1 ~ 



2k 0 



(ii) 



or 



k 0 = - + (x 1 - x 0 (Oo)) (x 2 - x 0 (9 0 )) + (Vi - Vo (Oo)) (1/2 - 2/o (Oo)) 



by (10). 

Solving for ko, and substituting this value into Claim 4, 5, or 6, depending on circumstances, determines 
the position of Target T. In fact, 



t = ti - j = h - 1 (/ici/l/co + t/ic) 2 ) 
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6. Calculation Algorithm 

To summarize, given SDU positions (xi , j/i), (x 2 , 2/2), and (X3, 2/3), and observed data {£1 , t 2 > £3}, fc he following 
algorithms give the Calculated Target Position (xq,2/o) of Target T. 



7. Common Elements 

A. Calculate f\ = h - £ 2 and f 2 = t\ - £3. 

#2 2 / 2 ' 

B. Calculate di = y (xi - x 2 ) + (1/1 - y 2 ) and d 2 = y (zi - z 3 ) + (2/1 - 2/3) 



8. When y x ^ y 2 , 2/1 # 2/3, and d< > c, for i = 1,2 

1. Calculate m x = - f^^V 

2. Calculate m 2 = - (^^). 

o tt . _ /ic A /4«+(/ l c) 2 +(»g-x;) + (yg-y?) 
USe 01 " 2(» a - Vl ) ' 

4 - USe > 2 = — * S ' 

5. Solve 

= g - (AC) 2 f (*1 - (dS=fc)) (*3 " (dfefc)) 

2 1 + (* - ("1 (dS=fc) + *)) (» - (m, (Jfe=fe) + 6,)) 

for k; call this value kq. 

6. Evaluate 61 and 6 2 with /c = /co; call these values /?i and /? 2) respectively. 

7. Then 

(2:0,2/0) = ( — — — ,mi [ — — —J + A J 
\mi - m 2 \mi -m 2 / / 



9. When 2/1 = 2/2 , 2/i 7^ 2/3 > and d< > c, for i = 1, 2 



1. Use ri = ^ C y^H7^|H-x?) 



2 - Use r2 = 2(5^ 



/ 2 c ((A +h)e± 4k + (f 1 c) 2 \ + (x 2 - x?) + (j/3 2 - y 2 ) 



3. Solve 



cP — ( f c)^ 

K = — ^ + ( Xl " ( X2 ~ + ^ ~ T ^ ^ 2 " ^ 



for k\ call this value kq. 

4. Evaluate r x and r 2 with ac = /cq; call these values 71 and 72, respectively. 

5. Then 

teo,yo) = (71,72) 



I 
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10. When yi ^ y 2 , 2/1 = 2/3> and di > |/i|c, for i = 1,2 
1 Use = /2c 2 (/i+/ 2 )±>/4ko+(/ 1 c) 2 +(x|-x;) 

1 2(13— Xl) 



_ t ±/iV4«o + ( / lC 7+(^- g?) + {yi - y?) 

2 " ^ ^ = ^ ( - [he + / 2 )c± v^T^?) + (*§ - x?)) (-ft) ) • 

3. Solve 

(P — ( f c) 2 

« = — — 2^ — + ( Xl ~ n ) ( X2 " ri ) + ( yi " r2 ^ ~ r2 ) 
for k; call this value kq. 

4. Evaluate ri and r 2 with /c = kq\ call these values 71 and 72, respectively. 

5. Then 

(20,2/0) = (71,72) 
11. Calculated Target Position Examples 

Suppose (21,2/1) = (2,3), (22,2/2) = (4,2), and (23,2/3) = (3,1), and observed data t\ = 1.6, t 2 = 2, 
£ 3 = 1.75. For demonstration purposes only, let c = 1 (see Figure ). Then y x ^ y 2 , yi ^ 2/3, and di > \fi\c, 
for i = 1,2. 

1. We have /1 = *i - t 2 = 1.6 - 2 = -0.40, and / 2 = - * 3 = 16 - 1.75 = -0.15. 

2. Then d x = ^(2 - 4) 2 + (3 - 2) 2 = \/5 and d 2 = ^/(2 — 3) 2 H- (3 — l) 2 = VS. 

3. And mi = - (§5§) = 2. 

4. Then m 2 = - (f^f) = \. 

5. Also b\ = °- 40 >/ 4lc+ (~ 0 -^jK 4 — 2 ) + ( 2 — ll = ^V25k -f 1 — |. Note that the positive sign has been 
chosen. 

-0.15 ( (-0.40-0. 15) + v^K+(-0.40) 2 ) + (3 2 -2 2 )+(l 2 -3 2 ) 1lft - „ . m 

6. And 6 2 = ^ 2(1-3) — = MM + 5§o^ 25 * + 1 - Note that the P osi ' 

tive sign has been chosen. 

7. Solving for gives us 



26454174801 561 * n „ nnA * nn « , ™™,o. 

k 0 = , nnnnnn „ dnn + ^ Of ^ rn V292162046321 = 1.5369135 • • ■ 
18790926400 2348865800 



8. Therefore 



2 /27205811857 561 /— — — — 7-; 7 
/3i = — \/— — — — + oo/ ^„ V292162046321- - = -2.9976991 • • 
^ 25 V 751637056 93954632 2 



and 



1167 3 /27205811857 561 , 

B 2 = — — + — \/ + V292162046321 = 0.8235564 • • • 

^ 2 1600 200 V 751637056 93954632 



9. So 



{x0iVo) = ( AzA ?mi ( AzA.) = | (2.5475037--- ,2.0973083 ■■■) 

\rrii — m 2 \mi — m 2 / / 

to 7 decimal places (see Figure 8). 

Now suppose (21,2/1) = (2,3), (22,2/2) = (4,3), and (23,2/3) = (3, 1), and the observed data is t\ = 1.6, 
t 2 =2, £3 = 1.75. Note that 2/1 = y 2 and 2/1 7^ 2/3 > an d d» > |/<|c, for i = 1,2. For demonstration purposes 
only, let c = 1. 

1. We have f x = *i - 1 2 - 1.6 - 2 = -0.40, and / 2 = *i - t 3 = 1.6 - 1.75 = -0.15. 

2. Then di = ^(2 - 4) 2 + (3 - 3) 2 = 2 and d 2 = \/(2 - 3) 2 + (3 - l) 2 = 

3. This means r x = -° a W 4k+ ^~°^ +( 16 ~ 4 ) _ _l_ y j2§K -f 1 + 3. Note that the positive sign has been 
chosen. 
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* T 

\ 

\ ' 



III 



Figure 8: Example Target Calculation Result when j/i ^ 3/2, V\ ^ J/3, and di > \ fi\c, for i = 1,2. 



4. And 



7*2 = 



2(1=31 



0.15 ((-0.40 - 0.15) ± ^/4k + (-0.40) 2 ) + (9 - 4) + (1 - 9) 



^ - (-0.40^4* + (-0.40) 2 + (16 - 4)) (§E§) 

= \ (-0.15 * 0 .55 - 0.0 5x/4k: + 0.16 + 9) 

- 3567 _ _1 /O^k + 1 
~~ 1600 200 v^*^ 1 

Note that the positive sign has been chosen. 
5. Solving for kq gives us 

5958301793 1233%/23857734 

«o = 



7 



6. Therefore 



3769960000 



117811250 



71 = "777 



6109100193 1233V23857734 



25 V 150798400 



+ 



4712450 



+ 3 = 2.7414202 



and 



72 = 



3567 



1 



1600 200 V 150798400 



6109100193 1233^23857734 n 

+ 1 = 2.1975504 



4712450 



7. So 



(zo.ito) = 1(2.7414202 2.1975504- »Q 
to 7 decimal places (see Figure 9). 



12. Algorithm Implementation Summary 

( TPR Calculation 16 ) 

(TPR Example #1 SDU Coordinates 18) 

(TPR Example #2 SDU Coordinates 23) 

( ELE Calculation 34 ) 

(ELE Example #1 SDU Coordinates 36) 

(ELE Example #2 SDU Coordinates 39) 

( Mitigations For Reflections 50 ) 

( Coverage Area Optimization 58 ) 



I 
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A 
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•III 



Figure 9: Example Target Calculation Result when yi = ?/2> Vi i 1 2/3, and d{ > \ fi\ c, for i = 1,2. 



13. Target Position Report Calculation 

A Target Position Report is the calculated (x, y) position value based on the arrival times detected at the 
SDU's, and the position information for the SDU's involved in the calculation. There are no adjustments 
made for the uncertainty in the arrival time data, as such issues are addressed by the Error Likelihood 
Ellipse methods. 

The principal calculation policy in the PQIC OSMMTS System is to base a Target Position Report on 
the choice of SDU arrival times that maximizes the likelihood that the calculated Target Position Report 
is as accurate, i.e., close to the actual target position, as possible. Since the PQIC OSMMTS SDU's do 
not change position as a function of time when the system is operating, this optimization may be achieved 
by implementing the PQIC OSMMTS Target Position Report Standard Methodology, which means pre- 
establishing the calculation algorithm that the system will utilize when a particular combination of SDU 
arrival times are used to calculate the Target Position Report. 

When error conditions are calculated, e.g., in the Error Likelihood Ellipse, the same algorithm set shall 
be used with appropriately adjusted arrival time data. 



14. Target Position Report Standard Methodology 

The PQIC OSMMTS Target Position Report Standard Methodology for calculating a Target Position Report 
begins with the use of a position-specific combination of the SDU's taken three at a time. Based on this 
position information, the Target Position Report algorithm calculates a (x,y) position value for the given 
hjhih values. This algorithm is specific to the combination of SDU's that receive the arrival time data. 
Therefore, for n-many SDU positions, there must be (3) -many algorithms pre-established for use by the 

Target Position Report calculation device. For example, if there are 5 SDU positions, then ( 3=10 ) separate 
algorithms must be established to calculate a Target Position Report regardless of which three SDU's are 
involved in the calculation. 

The set of algorithms generated for a particular set of SDU positions differ only by the number of terms 
and the coefficients of the terms involved in the sum. An effective organizing system that addresses the 
numerical issues related to such a complicated sum may be found later in this document. 
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15. Target Position Report MAPLE Algorithm 

Given the position information of the SDU's that receive the arrival time data, the following MAPLE algo- 
rithm calculates the Target Position Report based on the arrival times. For the purposes of this document, 
only example algorithms shall be given based on example position information. The completely generic 
form of the Target Position Report MAPLE Algorithm shall be documented in a separate PQIC OSMMTS 
technical memorandum. 

In general, each term of the Target Position Report calculation has the form 

(sign) x (coefficient) x (ratio) x ^ x t^ 1 x t 2 P2 x * 3 Ps 

for pi = 0, 1,2, where po = p\ + P2 + P3. The sign is either positive or negative 1, the coefficient is a 
positive integer, and the ratio is either 1 or a ration of previously calculated subexpressions. If the coefficient 
is allowed to be nonnegative, i.e., take on the value 0, then, in general, each Target Position Report may be 
expressed as the sum of the same number of such terms, where the number of such terms depends on the 
SDU position information. 

The Target Position Report terms may be organized into convenient subassemblies for calculation 
efficiency by using the tables found below as worksheets, as demonstrated later in this document. 

16. (TPR Calculation 16) = 
osmmts: = proc(nl , n2, n3) 

local fl , f2 , dl , ml , m2 , bl , b2 , A;, kOx , betal , beta2 ; 

globalc, xO , yO , xl , yl , x2 , y2 , x3 , y3 , si , s2 , s3 , kO , cthO , thO ; 

options i Copyright 2003PQ I Consulting A 11 Rights Reserved c ; 

descnp^on"OSMMTS u Target u Position u Report u aiid u Supporting u Inf ormation"; 

fl : = nl - n2\ f2: = nl - n3\ 

dl : = sqrt((x2 - xl ) © 2 + (y2 - yl ) 0 2); 

ml: = -((x2 -xl)/(y2 - yl))] 

m2: = -({x3 -xl)/{y3 -yl)); 

bl : = {fl * c * sqrt(4 * k + (fl * c) © 2) +{x2 © 2 - xl © 2) + (y2 © 2 - yl © 2))/(2 * (y2 - yl )); - 
b2: = (f2*c*((fl +f2)*c+sqrt(4*k+(fl *c)©2)) +{x3 ©2-xl ©2) + (y3 ©2- yl ©2))/(2* (y3 - yl )); 
kOx: 

{solve {k = ({dl @2-(fl *c)©2)/2) +(xl - ((b2 -bl)/(ml - m2))) *(x2 - ((b2 -bl)/(ml - m2))) 
+(yl - (ml * ((b2 - bl )/(ml - m2)) + bl )) *(y2 - (ml * ((b2 - bl )/(ml - m2)) + bl )), k)° }; 

if (nops(k0x) = 0) then 

kO: = undefined ; 

elif (nops(kOx) > 1) tfien 

if (&0x[2] — conjugate (A;0x[l])) tfien 

A;0: = undefined] 

else 

A;0: = mm(op(l, fcftr), op(2, fcOx)); 

end if ; 

else 

= Jbte[l]; 
end if ; 

ieiai: = eval(bl ,k = k0)\ 
beta2: = eva/(62,A; = fc0); 
xO: = (6e£a2 - betal)/(ml - m2)\ 
yO : = ml * xO + 6e£a! ; 

c£A0: = 1 - ((dl ®2-(}l * c) © 2)/(2 * k0))\ 
thO: = arccos(cthO); 

si : = sqrt((xl - xO) © 2 + (yl - t/0) © 2); 
52: = 5grt((x2 -x0)©2+(j/2 -y0)©2); 
55: = sqrt((x3 - xO) © 2 + (y5 - y0) © 2); 
return x0, y0, &0, &0x, nop5(A;0x), c£A0, th0,sl , s2>s3\ 
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endproc; 

This code is used in chunk 12. 



17. Target Position Report Example #1 

18. (TPR Example #1 SDU Coordinates 18) = 
(xl,yl): = (2,3); (x2,y2): = (4,2); (x8,yS): = (3,1); 

This code is used in chunk 12. 



19. Target Position Report MAPLE Calculation For Example #1 

> c:=l; (xl,yl) :=(2,3) ; (x2,y2) :=(4,2) ; (x3,y3) :=(3,1) ; 

c:= 1 
xl, yl := 2, 3 

x2, y2 := 4, 2 

atf, y3 := 3, 1 

> osmmts(tl,t2,t3) [1] ; 

20. X- Value Of Target Position Report For Example #1 

--(ft/ - ft?) c((2 « - t2 - tS) c + ((-(50 - L ~ 26c 4 tl 3 t2 - 30 tl c 4 ftS 3 
+ 16c 4 tl 2 t2 2 - 30 tl c 2 ft? - 12c 4 ft2 3 ftf + 26c 2 tl t2 + 70c 4 ftf 2 ft? 2 

+ 10 (5 a c< *' + 10 « c* <2 - 8 O c« If - %i a *f ts + "*<»°'» 2 

4%4ft9 c 3 ftg 2 4 4 ° 5%4cfc? %4 c ft/ 5 %4 c 3 ft? 3 4%4cft5 

+ %i +5c w ""^i - + _ %i %r~ + — %r~ 

2%4c 3 ft/ 3 

+ 57- 18c 2 ftg 2 + 2c 2 tl 2 + 46c 4 ft/ 2 /? ftg + 4c 4 tl t2 2 t3 

70 1 

- 50c 4 tl tS t3 2 + 10c 2 t3 2 + 22c 4 ft/ 4 + 5 c 4 ft§ 4 )/(%l) - (50 - 26c 4 tl 3 t2 

- 30 tl c 4 ft? 3 + 16c 4 tl 2 t2 2 - 30 ft/ c 2 ft? - 12 c 4 t2 3 tl + 26 c 2 /i ft2 + 70c 4 ft/ 2 ft? 2 

+ 10 t3 2 c 4 ftg 2 + 10 13 c 2 ft? - 8 13 c 4 t2 3 - 62 c 4 tl 3 t3 + 10 ft? 3 c 4 ft? 
4%3c 3 ft/ t2 2 7 %3 c 3 t J 2 fcg 13 %3 c 3 » 2 ft? 6 %3 c 3 t J fcg 1 3 

%1 + %1 %1 %1 

, 16%3ft/c 3 ft? 2 4%3«5c 3 i2 2 %3ft? 2 c 3 ft2 „ A . %3ctl 
+ ^ + ^ t^z + 5c 4 ft? 4 + 



%1 %1 %1 %1 

5%3c 3 ft? 3 4%3cftS 2%3c 3 ft/ 3 5%3cft? 



+ 46c 4 tl 2 t3 t2+Ac A tl t2 2 t3 - 50c 4 tl t2 t3 2 + 10c 2 ft? 2 + 22c 4 ft/ 4 + 5c 4 ft? 4 )/( 

%1)) + c 2 tl 2 - 2c 2 tl t2 + c 2 ft2 2 )( x /2)) + E + L tl _ t2 ) c((-(50 

6 3 

- 13%4 ^/ J - - 26c 4 tl 3 t2 - 30 tl c 4 i3 3 + 16c 4 £1 2 t2 2 - 30 W c 2 IS 

- 12 c 4 t2 3 tl + 26 c 2 « £B + 70 c 4 *1 2 2 + 10 tS 2 c 4 £2 2 + 10 tS c 2 t2 - 8 *3 c 4 i2 3 
%4 2 c 3 fcg 16 %4 « c 3 2 4 %4 c 3 « *2 2 7 %4 c 3 *1 2 12 _ 4 , n 

5gi — + %i %^- + — -%i 62cUl 13 

- 6 %4 ^ tf g + 10 tS 3 cH2 + A %4 ^ ^ 2 + 5 c 4 ^ 4 - 
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%Actl 5%4c 3 t3 3 4%4ctg 2%4c 3 tl 3 ,„ - o2 „ 2 <2 

+ %T- + -%T- + ^I w^tf' + a^" 1 

+ 46c 4 tl 2 t3 t2 + 4c 4 tl t2 2 t3 - 50c 4 tl t2 t3 2 + 10c 2 t3 2 + 22c 4 tl 4 + 5c 4 fg 4 )/( 
%1) - (50 - 26c 4 tl 3 t2 - 30 tl c 4 f$ 3 + 16c 4 tl 2 t2 2 - 30 tl c 2 t$ - 12c 4 tg 3 tl 

+ 26c 2 tl {2 + 70c 4 tl 2 t3 2 + 10 13 2 c 4 tg 2 + 10 13<? t2 - 8 13 c 4 fg 3 - 62c 4 tl 3 13 

lin , o34j „ 4%3c 3 tltg 2 7%3c 3 tl 2 fg U%3c 3 tl 2 t3 
4- 10 £j c £<; — H 

%i %i %i 

6 %3 c 3 tl tg t.9 16 %3 tl c 3 t3 2 4 %3 tg c 3 tg 2 %3 t3 2 c 3 tg 4 4 

%i + %i + %i %i +5c tf 

%3ctl 5%3c 3 t.? 3 4%3ctg 2_%3c 3 J£^ _ 5%3ct5 _ 2 2 
+ %1 %1 + %1 + %1 %1 18 C * 2 

+ 2c 2 tl 2 + 46c 4 tl 2 t3 t2 + 4c 4 tl t2 2 t3 - 50c 4 tl t2 t3 2 + 10c 2 t3 2 + 22c 4 tl 4 
+ 5c 4 tg 4 )/(%l)) + c 2 tl 2 -2c 2 tlt2 + c 2 tg 2 )* 1 / 2 ) 
%1 := 5c 2 tg 2 + 2c 2 tl 2 + 5c 2 t3 2 - 9 - 2 tl c 2 t3 - 2c 2 tl t2 - 8t3c 2 t2 

%2 := 50 - c 6 tl 2 t2 4 + 12c 6 tl 5 t2 + 6c 6 tl 3 tg 3 - 13c 6 tl 4 tg 2 + 12c 6 tl* t3 - c 6 tl 
+ 6c 6 ti 3 t3 3 - 13c 6 ti 4 tS 2 - 16c 4 tl 3 t2 - 14c 6 tl 2 t3 3 t2 + 34c 6 tl 3 t3 t2 2 

- 34c 6 tl 4 t3 t2 - t3 2 c 6 t2 4 - 2 t$ 3 c 6 fg 3 - tS 4 c 6 t2 2 + 34c 6 tl 3 tS 2 fg 

- 30 tl c 4 t3 3 + 8 c 4 tl 2 fg 2 - 30 c 6 tl 2 f3 2 fg 2 - 20 tl c 2 t$ - 10 c 4 tg 3 tl 

- 14c 6 tl 2 t2 3 tS + 20c 2 tl t2 + 10c 6 tl t2 2 t3 3 + 2c 6 tl t2 4 t3 + 10c 6 tl t2 3 t3 2 
+ 68 c 4 tl 2 t3 2 + 2 tl c 6 iS 4 t2 + 8 2 c 4 fg 2 + 50 t3 c 2 t2 - 10 f3 c 4 t2 3 - 4 c 6 f 1 6 

- 56 c 4 tl 3 t3 + 10 t3 3 c 4 fg + 5 c 4 f3 4 - 35 c 2 fg 2 + 32 c 4 tl 2 t3 t2 + 14 c 4 tl t2 2 t3 

- 46c 4 tl t2 t3 2 - 15c 2 tS 2 + 18c 4 tl 4 + 5c 4 tg 4 

%3 := 26c 3 tl 2 t3 - 14c 3 tl 2 tg + 8c 3 tl t2 2 - 32 tl c 3 t3 2 + 2 t3 2 c 3 t2 - 8 tSc 3 tg 2 
-8cfg + 10c 3 t3 3 -2cfl + 10cf3 + 12c 3 fl tg f3 - 4c 3 tl 3 - 6\/9o2 

%4 := 26c 3 tl 2 t3- 14c 3 tl 2 tg + 8c 3 tl t2 2 - 32 tl c 3 t3 2 + 2 t3 2 c 3 t2-8 13 c 3 tg 2 

- 8c tg + 10c 3 t3 3 - 2c tl + 10cf3 + 12 c 3 tl t2 t3 - 4c 3 tl 3 + 6\/%2 

21. y- Value Of Target Position Report For Example #1 

> osmmts(tl,t2,t3) [2] ; 

1 13%4r 3 tl 2 i? 

--(tl - t3) c((2 tl-t2- t3) c + ((-(50 - „" - 26c 4 tl 3 t2 - 30 tl c 4 t 
o Vol 

+ 16c 4 tl 2 fg 2 - 30 tl c 2 t3 - 12c 4 tg 3 tl + 26c 2 tl tg + 70c 4 tl 2 t3 2 

+ 10t 3 2 c 4 tg 2 + 10 t*c 2 tg -8t3c 4 tg 3 - %4 * + 16 

_ 4%4cHlt2 2 + 7%4c 3 U 2 t2 _ 62 c4 ^ 3 t3 _° 6%4^tttf' + 1Q ^ fi4 tf 

4%4t5c 3 tg 2 4 4 5%4ct<? %4 c tl 5 %4 c 3 t.? 3 4 %4 c tg 
+ %i +5c u ~~%T~ + %1 +_ %T~ 

+ 2 1°, * 1 - 18c 2 tg 2 + 2c 2 tl 2 + 46c 4 tl 2 f3 tg + 4c 4 tl fg 2 f3 

70 1 
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- 50c 4 tl t2 t3 2 + 10c 2 fg 2 + 22c 4 tl 4 + 5c 4 tg 4 )/(%l), -(50 - 26c 4 tl 3 t2 

- 30 tl c 4 t3 3 + 16c 4 tl 2 t£ 2 - 30 tl c 2 t3 - 12 c 4 t2 3 tl + 26c 2 tl t2 + 70 c 4 tl 2 fg 2 
+ 10 fg 2 c 4 fg 2 + 10 fg c 2 fg - 8 13 c 4 t2 3 - 62 c 4 tl 3 fg + 10 13 3 c 4 fg 

4 %3 c 3 tl fg 2 7 %3 c 3 t J 2 fg 13 %3 c 3 tl 2 t3 6 %3 c 3 tl t2 t3 
%l + %1 %1 %1 

16 %3 tl c 3 f g 2 4 %3 fg c 3 fg 2 %3 fg 2 c 3 t2 4 4 %3ctl 
+ %1 + %1 %1 +5 ° M + %1 

5%3c 3 fg 3 4%3cfg 2%3c 3 tl 3 5%3cfg , a , „, „ - <2 

%i + ^- + — %i %i 18c2 ^+ 2c2 " 2 

+ 46c 4 tl 2 fg t2 + 4c 4 tl fg 2 fg - 50c 4 tl fg t3 2 + 10c 2 fg 2 + 22c 4 tl 4 + 5c 4 fg 4 )/( 

961)) + c 2 fl 2 - 2 c 2 tl t2 + c 2 fg 2 )* 1 / 2 )) + ^ + \{tl - t2) c((-(50 

6 6 

- lS%i ^J l2tS _ 26c 4 tl 3 t2 - 30 tl c 4 t3 3 + 16c 4 tl 2 t2 2 - 30 tl c 2 tf 

- 12 c 4 fg 3 tl + 26 c 2 tl tg + 70 c 4 tl 2 fg 2 + 10 13 2 c 4 fg 2 + 10 13 c 2 t2 - 8 fg c 4 fg 3 
%4 fg 2 c 3 fg 16 %4 ^ic 3 fg 2 4 %4 c 3 tl t2 2 7 %4 c 3 f 1 2 tl 4 3 

%i + %i %i + %i 

6%4c 3 flfgfg in . ?3c 4,g , 4%4fgc 3 fg 2 . 5%4ctg 

%i +1UW ° w+ %1 " %1 

%4ctl 5%4c 3 fg 3 4%4cfg 2%4c 3 tl 3 1Q - n2 „ 2i<2 

+ %r- + -^r- + — %i — ^ 2 + 2c 2 tl 2 

+ 46c 4 tl 2 t3 t2 + 4c 4 tl t2 2 t3 - 50c 4 tl t2 t3 2 + 10c 2 tg 2 + 22c 4 tl 4 + 5c 4 fg 4 )/( 
%1), -(50 - 26c 4 fl 3 t2 - 30 fl c 4 t3 3 + 16c 4 fl 2 tg 2 - 30 tl c 2 tg - 12c 4 t2 3 tl 

+ 26 c 2 tl fg + 70 c 4 f 1 2 tg 2 + 10 tg 2 c 4 fg 2 + 10 tg c 2 fg - 8 tg c 4 t2 3 - 62 c 4 tl 3 tg 

3 4 4%3c 3 tl tg 2 7%3c 3 tl 2 t2 13 %3 c 3 tl 2 tg 
+ 10 tg c fg %1 + %1 %1 

6 %3 c 3 tl tg tg 16 %3 tl c 3 tg 2 4%3fgc 3 fg 2 %3 tg 2 c 3 tg 4 4 
%1 + %1 + %1 %1 +5c " 

%3ctl 5%3c 3 tg 3 4%3cfg 2%3c 3 Jl 3 _ _ 5%3ctg _ 2 2 
+ %1 %1 + %1 + %1 %1 

+ 2c 2 tl 2 + 46c 4 tl 2 tg tg + 4c 4 tl tg 2 tg - 50c 4 tl tg tg 2 + 10c 2 tg 2 + 22c 4 tl 4 
+ 5c 4 fg 4 )/(%l)) + c 2 fl 2 - 2c 2 tl fg + c 2 fg 2 )( 1 / 2 ) 

%1 := 5c 2 tg 2 + 2c 2 tl 2 + 5c 2 tg 2 - 9 - 2 tl c 2 tg - 2c 2 tl tg - 8 fgc 2 tg 

%2 := 50 - c 6 tl 2 tg 4 + 12c 6 tl 5 tg + 6c 6 tl 3 tg 3 - 13c 6 tl 4 tg 2 + 12c 6 tl 5 tg - c 6 tl 2 tg 4 

+ 6c 6 tl 3 tg 3 - 13c 6 tl 4 tg 2 - 16c 4 tl 3 t2- 14c 6 tl 2 tg 3 tg + 34c 6 tl 3 tg tg 2 

- 34 c 6 tl 4 t3t2- tg 2 c 6 tg 4 - 2 tg 3 c 6 tg 3 - tg 4 c 6 tg 2 + 34 c 6 tl 3 tg 2 tg 

- 30 tl c 4 tg 3 + 8 c 4 tl 2 tg 2 - 30 c 6 tl 2 tg 2 tg 2 - 20 tl c 2 tg - 10 c 4 tg 3 tl 

- 14 c 6 tl 2 fg 3 tg + 20 c 2 tl tg + 10 c 6 tl tg 2 fg 3 + 2 c 6 tl fg 4 tg + 10 c 6 tl tg 3 tg 2 
+ 68c 4 tl 2 tg 2 + 2 tl c 6 tg 4 tg + 8 fg 2 c 4 tg 2 + 50 tg c 2 tg - 10 fgc 4 tg 3 - 4c 6 tl 6 

- 56c 4 tl 3 tg + 10 fg 3 c 4 fg + 5c 4 tg 4 - 35c 2 tg 2 + 32c 4 tl 2 tg tg + 14c 4 tl tg 2 tg 

- 46c 4 tl tg fg 2 - 15c 2 tg 2 + 18c 4 tl 4 + 5c 4 tg 4 



PQIC Documentation — OSMMTS Patent Application Technical Documentation 



%3 := 26c 3 tl 2 t3 - 14c 3 tl 2 t2 + 8c 3 tl t2 2 - 32 tl c 3 t3 2 + 2 2 c 3 t2 - 8 c 3 t2 2 
-8c*2 + 10c 3 t3 3 -2ctl + 10ctf + 12c 3 *2 *3-4c 3 <1 3 -6\/%2 
%4 := 26c 3 « a - 14c 3 tl 2 12 + 8c 3 *i fcg 2 - 32 *i c 3 t3 2 + 2 *3 2 c 3 t2 - 8 *3 c 3 t2 2 
-8c ^4- 10c 3 tf 3 -2ci/ + 10 ct3 + 12c 3 *i *2 tS - 4c 3 *1 3 + 6 y/%2 

22. Target Position Report Example #2 

23. (TPR Example #2 SDU Coordinates 23) = 
(xl,yl): = (3,2); (x2,y2): = (2,4); (atf.yJ): = (1,3); 

This code is used in chunk 12. 



24. Target Position Report MAPLE Calculation For Example #2 

> c:=l; (Xl,yl):=(3,2);(x2,y2):=(2,4);(x3,y3):=(l,3); 

c:= 1 
xl, yl := 3, 2 

x2, t/2 := 2, 4 

x3, y3 := 1, 3 

> osmmts(tl,t2,t3) [1] ; 

25. X- Value Of Target Position Report For Example #2 

~(tl - t3) c((2 tl -t2- t3) c + (4min(-(50 + 26 c 2 tl t2 - & t3 c 4 i£ 3 - 30 tl c 2 t3 

+ 10 tS c 2 t2 + 10 t3 3 c 4 *g - 6 %3 C ^ * * + 4 %3 1^ *' - %3 %f 12 

7%3c 3 tl 2 t2 4%3c 3 tlt2 2 16 %3 ti c 3 *3 2 13 %3 c 3 » 2 t3 
+ %1 %1 + %1 %1 

- + io t 5 2 c 4 t2 2 - 30 « c 4 t3 3 + 70c 4 ii 2 tf 2 - 12c 4 t2 3 tl 
vol 

+ 16c 4 tl 2 t2 2 + 5c 4 <3 4 - 18c 2 *2 2 + 2c 2 t/ 2 +46c 4 tJ 2 t2 + 4c 4 *2 2 t3 

- 50c 4 «J fcg iS 2 + 10c 2 t3 2 + 22c 4 tl 4 + 5c* t2 4 )/{4%l), -(50 + 26c 2 tl t2 

4 %4 f$ c 3 *2 2 



- 8 13 c 4 i2 3 - 30 tl c 2 (3 + 10 t3c 2 t2 + 10 13 3 c 4 «2 + 



%1 



_%U3^J2 7%4c 3 tl 2 t2 4%4<? tlt2 2 16 %4 tl c 3 13 2 

%1 + %1 %1 + %1 

13%4c 3 f/ 2 t3 6%4c 3 Wt2t3 „„ . ~ „ „„ . - „ „, , „, 

57- 62 c 4 tl 3 t3 - 26c 4 tl 3 t2 + 10 W 2 c 4 fc2 2 

/ol %1 

5%4c 3 i3 3 4%4ct2 5%4ct3 %4ctl 2%4c 3 tl 3 _ „ „, 

-— %r- + %r- + ^i- + ^i — 30tlcts 

+ 70c 4 tl 2 t3 2 - 12c 4 t2 3 tl + 16c 4 *i 2 t2 2 + 5c 4 tf 4 - 18c 2 t2 2 -Vic 2 tl 2 

+ 46c 4 W 2 t3 t2 + 4c 4 tl t2 2 t3 - 50c 4 tl t2 t3 2 + 10c 2 t3 2 + 22c 4 tl 4 +5c 4 t2 4 )/( 

4 961)) + c 2 2 - 2 c 2 t2 + c 2 fS 2 )* 1 / 2 *) + ^ + Utl - *2)c(4min(-(50 

6 6 
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+ 26 c 2 tl t2-St3c 4 t2 3 - 30 tl c 2 tf + lOtSc 2 tf + 10« 3 c 4 « 

6 %3 c 3 tl tg tff 4 %3 ft? c 3 ^ 2 963 2 c 3 tg 7%3c 3 tl 2 t2 

961 * %1 961 + 961 
4%3c 3 tlt2 2 16%3tlc 3 t3 2 U%3c 3 tl 2 t3 ^ 4 ^ 3i0 963 c« 
%i + %i %i 62c 4 tl 3 t 3 + — 

2%3c 3 tl 3 _ 5%3c 3 t3 3 4 3 4963ctg 5%3ctf 
+ 961 961 2bC " ^ + %l 961 

+ 10 t3 2 c 4 t2 2 - 30 tl c 4 t3 3 + 70 c 4 tl 2 2 - 12 c 4 t2 3 tl + 16 c 4 tl 2 tg 2 + 5 c 4 t,? 4 

- 18c 2 t2 2 + 2c 2 tl 2 + 46c 4 tl 2 tf tg + 4c 4 tl tg 2 t3 - 50c 4 tl tg t3 2 + 10c 2 t3 2 

+ 22 c 4 tl 4 + 5 c 4 tg 4 )/(4 961), -(50 + 26 c 2 tl tg - 8 tS c 4 tg 3 - 30 tl c 2 tS 

,n « ,n.o3 4 4 964 W c 3 tg 2 %4 13 2 c 3 t2 7%4c 3 tl 2 t2 

+ 10 t3 c 2 tg + 10 t3 3 C 4 tg + 5^; 57; + 57- 

Tol 7ol /ol 

4%4c 3 tl t2 2 16 %4 tl c 3 t.3 2 13 %4 c 3 tl 2 tS 6 %4 c 3 tl tg ft? 

961 + 961 961 961 

_ , . n _ 4 o ^ „ 2 4 - 5%4c 3 t3 3 4%4ctg 5%4ct3 

- 62 c 4 tl 3 t3 - 26c 4 tl 3 tg + 10 t,? 2 c 4 tg 2 r- + — 5- — — 

%1 %1 %1 

%4 c tl 2 %4 c 3 tl 3 .0 ^00 j o ^00 

4- ° + 30 tl c 4 t3 3 + 70c 4 tl 2 tf 2 - 12c 4 tg 3 tl + 16c 4 tl 2 t2 2 

Vol Vol 

4- 5 c 4 t3 4 - 18c 2 tg 2 + 2 c 2 tl 2 + 46c 4 tl 2 tS tg + 4c 4 tl tg 2 t3 - 50c 4 tl fig t3 2 
+ 10c 2 tf 2 + 22 c 4 tl 4 + 5c 4 tg 4 )/(4%l)) + c 2 tl 2 - 2 c 2 tl tg + c 2 tg 2 )<V2) 

961 := 2c 2 tl 2 + 5c 2 t3 2 + 5c 2 tg 2 - 9 - 8 c 2 tg - 2 tl c 2 t3 - 2c 2 tl ft* 

962 := 50 + 20c 2 tl t2 - 10 tS c 4 tg 3 - 20 tl c 2 tS + 50 t3 c 2 tg + 10 13 3 c 4 tg - 56c 4 tl 3 tS 

- 16c 4 tl 3 12 + 8 t3 2 c 4 tg 2 - 30 tl c 4 t2 3 + 68c 4 tl 2 13 2 - 10c 4 tg 3 tl + 8c 4 tl 2 t2 2 

+ 5c 4 t3 4 - 35c 2 tg 2 + 32c 4 tl 2 tS t2 + 14c 4 tl tg 2 tf - 46c 4 tl t2 t3 2 - 15c 2 t3 2 
+ 18c 4 tl 4 + 5c 4 tg 4 - 4c 6 tl 6 - 34c 6 4 t2 + 34c 6 il 3 « t2 2 - 14c 6 ii 2 « 3 t2 

- 30c 6 il 2 « 2 t2 2 + 34c 6 tl 3 tS 2 fg - 14c 6 tl 2 ig 3 tS + 6c 6 ii 3 i3 3 + 12c 6 il 5 tf 

- c 6 il 2 « 4 - 13c 6 tl 4 t2 2 + 6c 6 tl 3 *g 3 + 12c 6 tl 5 tg - c 6 tl 2 tg 4 - t3 4 c 6 <g 2 

- 2 t3 3 c 6 ig 3 - t5 2 c 6 ig 4 - 13 c 6 tl 4 t3 2 + 10 c 6 tl tg 3 tf 2 + 2 c 6 tl t2 4 tS 
+ 10 c 6 tl t2 2 t3 3 + 2 tic 6 t3 4 t2 

963 := 26c 3 tl 2 - 14c 3 tl 2 tg + 8c 3 tl t2 2 - 32 tl c 3 t3 2 + 2 t3 2 c 3 tg - 8 t3c 3 tg 2 

- 4c 3 tl 3 + 10c 3 t3 3 + lOctf - 8ctg - 2ctl + 12c 3 tl tg t5 -f 6\/%2 

964 := 26 c 3 tl 2 t3 - 14 c 3 tl 2 t2 + 8 c 3 tl t2 2 - 32 tl c 3 t3 2 + 2 t3 2 c 3 tg - 8 15 c 3 tg 2 
-4c 3 tl 3 + 10c 3 t3 3 + 10ct5-8ctg-2ctl +12c 3 tl tg -6>/%2 

26. y- Value Of Target Position Report For Example #2 

> osmmts(tl,t2,t3) [2] ; 

-\{tl - tS) c((2 tl - tg - tf ) c + (4min(-(50 + 26 c 2 tl tg - 8 13 c 4 tg 3 - 30 tl c 2 t3 
6 

1Alfl 2Afl inxo3 4 ^ 6%3c 3 tlt2t3 t 4%3t3c 3 t2 2 %3t3 2 c 3 t2 
+ 10 15 c 2 tg + 10 13 3 c 4 t2 5- + 7— ^ 

70l 7ol 701 

7 963 c 3 tl 2 tg 4%3c 3 tl tg 2 16 963 tlc 3 t3 2 13 963 c 3 tl 2 t.9 
+ 961 961 + 961 %1 

62. 4 tl 3 tJ 1 %3c " l 2%3c ^ i3 5%3c3 ^ 26c 4 tl 3 tg l 4%3ct2 
b2ctltJ+ %1 + %1 2bc " ^ + 961 
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_ 5%3ctg + 1Q ^ 2 c4 tg2 _ 3Q ^ c4 ^ + ?o c4 tj2 ^ _ tg3 ^ 

+ 16c 4 tl 2 t2 2 + 5c 4 tg 4 - 18c 2 fg 2 + 2c 2 tl 2 + 46c 4 tl 2 tg t2 + 4c 4 tl tg 2 tg 
- 50c 4 tl tg t3 2 + 10c 2 tg 2 + 22c 4 tl 4 + 5c 4 *g 4 )/(4%l), -(50 + 26 <Ptl t2 

4 %4 tg c 3 t2 2 



-8t3c 4 t2 3 - 30 tl c 2 tg + 10 igc 2 tg + 10 tg 3 c 4 tg + 



%1 



_%4£ 2 _c 3 j£ 7%4c 3 tl 2 t2 4%4c 3 tlt2 2 16 %4 tl c 3 tg 2 

%1 + %1 %1 + %1 

13%4c 3 tl 2 fg 6%4c 3 tlt2t3 m . - „ M 4 <3 „ „, 4 „, 

^ rr- 62 c 4 tl 3 tg - 26c 4 tl 3 t2 + 10 fg 2 c 4 t2 2 

/ol %1 

5%4c 3 fg 3 4%4ct2 5%4c<5 %4cfl 2%4c 3 tl 3 

— %r— + ^i %^ + ^i" + — xi — 30tlcts 

+ 70c 4 tl 2 tg 2 - 12c 4 t2 3 tl + 16c 4 tl 2 t2 2 + 5c 4 t3 4 - 18c 2 fg 2 + 2c 2 tl 2 

+ 46c 4 tl 2 t3 t2+4c 4 tl t2 2 t3 - 50c 4 tl t2 t3 2 + 10c 2 t3 2 + 22c 4 tl 4 + 5c 4 tg 4 )/( 

4%1)) +c 2 tl 2 - 2c 2 tl t2 + c 2 tg 2 )^)) + H + i(tl - fg)c(4min(-(50 

6 3 

+ 26 c 2 tl t2 - 8 tg c 4 fg 3 - 30 tl c 2 tg + 10 13 c 2 t2 + 10 13 3 c 4 t2 
6 %3 c 3 tl tg t3 4 %3 tg c 3 tg 2 %3 fg 2 c 3 tg 7 %3 c 3 tl 2 t2 

%1 + %1 %1 + %1 

4%3c 3 tltg 2 16%3flc 3 tg 2 13 %3 c 3 tl 2 tg M 4 n %3ctl 

%1 + %1 %i 62c 4 tl 3 tg + -^- 

2%3c 3 tl 3 5%3c 3 tg 3 nc 4 <3 4%3ct£ 5%3ct£ 

+ %1 %1 6 ° " * 2+ %1 %1 

+ 10 fg 2 c 4 tg 2 -30 tlc 4 tf 3 + 70c 4 tl 2 tg 2 - 12c 4 tg 3 tl + 16c 4 tl 2 t2 2 + 5c 4 t3 4 
- 18c 2 tg 2 + 2c 2 tl 2 + 46c 4 tl 2 tg t2 + 4c 4 tl t2 2 t3 - 50c 4 tl tg tg 2 + 10c 2 fg 2 

+ 22c 4 tl 4 + 5c 4 tg 4 )/(4%l), -(50 + 26c 2 tlt2-8 tSc 4 tg 3 - 30 tl c 2 t3 

4 %4 13 c 3 tg 2 %4 tg 2 c 3 tg 7 %4 c 3 tl 2 tg 



+ 10 fg c 2 tg + 10 tg 3 c 4 tg + 



%1 %1 %1 



4 %4 c 3 tl tg 2 16 %4 tl c 3 tg 2 13 %4 c 3 tl 2 tg 6 %4 c 3 tl tg tg 
%1 + %1 %1 %1 

- 62c 4 tl 3 tg - 26c 4 tl 3 tg + 10 tg 2 c 4 t2 2 - + - 

%1 %1 /ol 

+ + 2% t C !" 3 " 30 tl c 4 tg 3 + 70c 4 tl 2 tg 2 - 12c 4 tg 3 tl + 16c 4 tl 2 tg 2 

Vol 701 

+ 5c 4 tg 4 - 18c 2 tg 2 + 2c 2 tl 2 + 46c 4 tl 2 tg tg + 4c 4 tl tg 2 tg - 50c 4 tl tg tg 2 
+ 10c 2 tg 2 + 22c 4 tl 4 + 5c 4 tg 4 )/(4%l)) + c 2 tl 2 - 2 c 2 tl tg + c 2 t2 2 )W 

%1 := 2c 2 tl 2 + 5c 2 tg 2 + 5c 2 tg 2 - 9 - 8tgc 2 tg - 2 tl c 2 tg - 2c 2 tl tg 

%2 := 50 + 20 c 2 tl tg - 10 fgc 4 tg 3 - 20 tl c 2 tg + 50 tgc 2 tg + 10 tg 3 c 4 tg - 56c 4 tl 

- 16c 4 tl 3 fg + 8 tg 2 c 4 tg 2 - 30 tl c 4 tg 3 + 68c 4 tl 2 tg 2 - 10c 4 tg 3 tl + 8c 4 tl 2 tg 2 

+ 5c 4 tg 4 - 35c 2 tg 2 + 32c 4 tl 2 fg tg + 14c 4 tl tg 2 tg - 46c 4 tl tg tS 2 - 15c 2 tg 2 
+ 18c 4 tl 4 + 5c 4 tg 4 - 4c 6 tl 6 - 34c 6 tl 4 tg tg + 34c 6 tl 3 tg tg 2 - 14c 6 tl 2 t3 3 t2 

- 30c 6 tl 2 tg 2 tg 2 + 34c 6 tl 3 tg 2 tg - 14c 6 tl 2 tg 3 tg + 6c 6 tl 3 tg 3 + 12c 6 tl 5 fg 

- c 6 tl 2 tg 4 - 13c 6 tl 4 tg 2 + 6c 6 tl 3 tg 3 + 12c 6 tl 5 tg - c 6 tl 2 tg 4 - tg 4 c 6 tg 2 

- 2 tg 3 c 6 tg 3 - tg 2 c 6 tg 4 - 13c 6 tl 4 tg 2 + 10c 6 tl tg 3 tg 2 + 2c 6 tl tg 4 tg 
+ 10 c 6 tl tg 2 tg 3 + 2 tl c 6 tg 4 tg 

%3 := 26c 3 tl 2 tg - 14c 3 tl 2 fg + 8c 3 tl tg 2 - 32 tl c 3 tg 2 + 2 tg 2 c 3 t2 - 8 tg c 3 tg 2 

-4c 3 tl 3 + 10c 3 tg 3 + 10ctg-8cfg-2cfl + 12c 3 fl fgfg + 6\/%2 

%4 := 26c 3 tl 2 tg - 14c 3 tl 2 tg + 8c 3 tl tg 2 - 32 tl c 3 t3 2 + 2 tg 2 c 3 tg - 8 tg c 3 fg 2 
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-4c 3 tl 3 + 10c 3 t3 3 + 10ct3-8ct2 - 2ctl + 12c 3 tl t2 t3 - 6a/%2 

27. Assembly Calculation Sequence For A Target Position Report 

The order of arithmetic steps to calculate a Target Position Report, given arrival times t\,t 2 ,h, does not 
depend on the coefficients found in a particular MAPLE calculation. In fact, a generic calculation sequence is 
available regardless of the values assigned to the (a*, For the purposes of this section only, the coefficient 
and pi values for the X- value of Target Position Report Example #1 shall be used as a demonstration of the 
calculation sequence that may be utilized to assemble a Target Position Report. This sequence is oriented 
to the anticipated MMIX assembly language implementation that forms the PQIC OSMMTS instantiation. 

%1 := 5c 2 t2 2 + 2c 2 tl 2 + 5c 2 tS 2 - 9 - 2 tl c 2 t3 - 2c 2 tl t2 - 8 t3c 2 t2 

%2 := 50 - c 6 tl 2 t2 4 + 12c 6 tl b t2 + 6c 6 tl 3 t2 3 - 13c 6 tl 4 t2 2 + 12c 6 tl 5 t3 - c 6 tl 2 t3 4 
+ 6c 6 tl 3 t3 3 - 13c 6 tl 4 t3 2 - 16c 4 tl 3 t2 - 14c 6 tl 2 t3 3 12 + 34c 6 tl 3 t3 12 2 

- 34 c 6 tl 4 t3t2 - t3 2 c 6 t2 4 - 2 tS 3 c Q t2 3 - t3 4 c 6 t2 2 + 34 c 6 *i 3 t3 2 t2 

- 30 tl c 4 t3 3 + 8c 4 *1 2 *2 2 - 30c 6 fi 2 fc? 2 fg 2 - 20 *i c 2 *5 - 10c 4 fcg 3 tl 

- 14 c 6 « 2 + 20 c 2 « fiB + 10 c 6 « *3 3 + 2 c 6 O *2 4 « + 10 c 6 U *2 3 t3 2 
+ 68c 4 « 2 t3 2 + 2 *1 c 6 4 *jS + 8 2 c 4 *2 2 + 50 c 2 £8 - 10 « c 4 *2 3 - 4c 6 £i 6 

- 56c 4 tl 3 t3 + 10 3 c 4 fcg + 5 c 4 *3 4 - 35c 2 t2 2 + 32c 4 il 2 *3 fcg + 14c 4 « t2 2 

- 46c 4 t2 t3 2 - 15c 2 2 + 18c 4 tl 4 + 5c 4 *2 4 

%3 := 26 c 3 « 2 - 14 c 3 tl 2 t2 + 8 c 3 *i - 32 tl c 3 2 + 2 t3 2 c 3 t2 - 8 *3 c 3 *2 2 
-8c*2 + 10c 3 t3 3 -2ctl +10c*3 + 12c 3 tl t2 t3 -4c 3 il 3 ~ 6 \/%2 



%4 := %3 + 12\/%2 



# 1 := _( 50 - 13 ^yj 1 2 13 - 26c 4 tl 3 tg - 30 « c 4 tf 3 

+ 16 c 4 tl 2 t2 2 - 30 tl c 2 t3 - 12 c 4 *2 3 M + 26 c 2 *i + 70 c 4 tl 2 t3 2 

+ io *» c 4 * » + io « c* * - s ts c* * > - %4t i 2c3tg + H %4 " c3 * 2 



%i 



%i 



4%4c 3 iJi2 2 7%4c 3 tl 2 t2 cn . - 6%4c i tlt2t3 ,„ , , „ 

+ — 62c 4 <i 3 tS — + 10 13 3 c 4 t2 



%1 



%1 



%1 
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4 %4 t3 c 3 tg 2 4 4 5%4ct3 %4ctl 5%4c? t3 3 4%4ctg 

+ %i +5c u + %i +— %r~ 

+ q ' c _ 18c 2 tg 2 + 2c 2 ti 2 + 46c 4 ti 2 t3 t2 +4c 4 ti tg 2 t5 
- 50c 4 ti tg t3 2 + 10c 2 tS 2 + 22c 4 ti 4 + 5c 4 tg 4 )/(%l) 



#2 := -(50- 26c 4 ti 3 tg 

- 30 ti c 4 W 3 + 16c 4 tl 2 t2 2 - 30 ti c 2 t3 - 12c 4 tg 3 ti + 26c 2 ti 42 + 70c 4 tl 2 t3 2 
+ 10 13 2 c 4 tg 2 + 10 13 c 2 t2-8 13 c 4 tg 3 - 62 c 4 ti 3 tS + 10 t3 3 c 4 tg 

4%3c 3 ti tg 2 7 %3 c 3 ti 2 12 13 %3 c 3 ti 2 6 %3 c 3 ti tg 
%1 + %1 %1 %1 

16 %3 ti c 3 t.9 2 4 %3 tff c 3 tg 2 %3 t3 2 c 3 t2 4 4 %3 c tl 

+ %! + %i %i +5c w + ~%T~ 

5%3c 3 f3 3 4%3cfg 2%3c 3 i/ 3 5%3cf3 ,„ 2 - „ 2 <2 
ST" + + %1 

+ 46c 4 2 tS tg + 4c 4 tl t2 2 t3 - 50c 4 tl t2 t3 2 + 10c 2 t3 2 + 22c 4 ti 4 + 5c 4 tg 4 )/( 
%1) 



#3 := c 2 tl 2 -2c 2 tlt2 + (? t2 2 



#0 := (2 tl -t2- tS) c 



#00 := -\{tl -t3)c 
6 



A := (#0 + V#l + #2 + #3) x #00 



#5:=|(« -tg)c 



13%4r 3 t/ 2 /!? 

#6 := -(50 - 1" - 26c 4 tl 3 tg - 30 ti c 4 t3 3 + 16c 4 ti 2 tg 2 - 30 ti c 2 t3 

vol 

- 12 c 4 tg 3 ti + 26 c 2 ti tg + 70 c 4 ti 2 t3 2 + 10 t3 2 c 4 tg 2 + 10 15 c 2 tg - 8 13 c 4 tg 3 
%4t3 2 <?t2 16%4tic 3 t5 2 4%4c 3 fi tg 2 7%4c 3 tl 2 12 



%1 %1 %1 T %1 



-62c 4 ti 3 t5 



6%4c 3 tlt2t3 1Q . 3 3 c 4, g , 4%4t3c 3 tg 2 4 4 5%4ctg 
%1 c W+ +5c W - 



I 
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, %4ctl 5%4c 3 tg 3 4%4ct2 2%4c 3 tl 3 , 2 , 2 

+ %! — + — %i + %i 18c ^ +2c " 

+ 46 c 4 ti 2 tS ti + 4 c 4 ti t2 2 t3 - 50 c 4 tl t2 t3 2 + 10 c 2 f3 2 + 22 c 4 ti 4 + 5 c 4 f2 4 )/( 
%1) 



#7 := -(50 - 26c 4 tl 3 t2 - 30 tl c 4 tS 3 + 16c 4 tl 2 t2 2 - 30 ti c 2 t3 - 12 c 4 t2 3 tl 
+ 26c 2 t2 + 70c 4 ti 2 t3 2 + 10 13 2 c 4 fg 2 + 10 t3c? t2 - 8 t$c 4 tS 3 - 62c 4 tl 3 t3 

4%Zc 3 tlt2 2 !%Z(?tl 2 t2 13%3c 3 ti 2 t3 

+ 10 « 3 c 4 tf 555: — + — 551 as 

6 %3 c 3 ti t2 t3 16 %3 ti c 3 t3 2 4 %3 13 c 3 f2 2 %3 t3 2 c 3 tg _ 4 o4 

%i + %i + %i %i +5c 45 

%3cti 5%3c 3 t3 3 4%3ct2 , 2%3(? tl 3 5%Zct3 2 2 

+ %i— + + — %I %i — 18c t2 

+ 2c 2 ti 2 + 46c 4 ti 2 f3 t£ + 4c 4 ti t2 2 t3 - 50c 4 ti t2 t3 2 + 10c 2 t5 2 + 22c 4 ti 4 
+ 5c 4 t2 4 )/(%l) 



#8 -c 2 ^ 2 -2c 2 ti t2 + c 2 t2 2 
B := (#5 x V#6 + #7 + #8) + #4 
So := vl + B 
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28. Error Distribution 

All data collected during process actions axe subject to error. The source or cause of the error varies by 
implementation, however, the OSMMTS, like all such processes, has several dominant sources of data error. 
There is mechanical error introduced by the physical limitations of the mechanical devices to precisely 
calculate the exact time of signal detection. There is propagation error from the natural corruption of the 
signal before detection. There is also rounding error due to the finite precision calculations utilized in the 
implementation, and algorithm error introduced by inefficient numerical methods. 

The PQIC OSMMTS may be used to specifically addresses rounding and algorithm errors by framing the 
implementation of all algorithms in assembly code, e.g., within a system that uses the MMIX 11 instruction set 
and memory model, even though this approach is not required in any particular OSMMTS implementation 12 . 
However, utilizing MMIX-based assembly code may be considered the OSMMTS "preferred embodiment," 
with regard to algorithm implementation, since all non-mechanical and non-propagation error issues are 
addressed by this approach. 

Finally, propagation error may be addressed, if not eliminated, by the optimal placement and mitigations 
for reflections methods documented in this memorandum. Mechanical errors may be addressed by the 
merchantability of the PCPU. 

The primary focus of the error distribution calculations is to account for the remaining effect of error 
in the arrival time data, regardless of source, after all precautions, optimizations, and mitigations have been 
used to minimize said error. 

Definition 8 The MLTarget Position is the calculated target position most likely to be correct given the 
arrival time data and the variability of the arrival times. It is found at the same time the Mo distribution 
is calculated (see 9). 

29. Derivation 

Suppose t\ and t 2 are not precisely measured, in the sense that 

t\ = t\ + E\ and <2 = *2 + £2 

where the t • are the observed values, and the U are the actual "true" times. Here, each e» is an error term. 
Suppose each 6i ~ N (0,a 2 ), for a common o > 0, and all such e» are IID. Then 

/• = *S-*5 

= (ti+ei)-(t 2 +e 2 ) 
= (*i - t 2 ) + (e x - £ 2 ) 

= / + (ci - £2) 

which means 

f*-f = ei-e 2 ~ AT (0,2a 2 ) 

Then 

5* — si = f*c and s\ - s 2 = fc 

mean 

W - *5) - («i - *) = (/•-/) c - N (0, 2cV) 

or 

Mi " «n) " (*2i - sn) = - (A* - fx) c ~ N (0, 2cV) 

(S22 ~ *ia) " (S22 - 52i ) = - US " /a) c " N (0, 2cV) 
Now from (1), we have 

d\ = s 2 u + s\x - 2sns 2 i cos 0 O 

and 

dy — S 1 ^ "h $21 """" ^^11*^21 COS 0Q 

11 MMIX (2009): A RISC Computer For The Third Millennium, as developed by Dr. Donald E.'Knuth of Stanford University, 
is documented at http://wv-cs-f acuity .Stanford. edu/"knuth/mmix. html. 

12 An OSMMTS implementation will work as well without MMIX as any other implemented OSMMTS as long as the error 
issues identified herein are addressed. 
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so 



s n + s 2i ~ 2 *o cos tf 0 = s{l + - 2k o cose o 



However, 



s ii + S 21 



2 \ \ 

U ^ci^Ko + ^ic) 2 )) +(| (-/ic±74«o + (/ic) 2 JJ 

i (4 (/ lC ) 2 + 8/co) 
(/ lC ) 2 + 2ko 



and 



*i 2 + «2 2 = (/rc) 2 + 2«s 

(/ lC ) 2 + 2k 0 (1 - cos0 o ) = (/i*c) 2 + (1 - cos^) 
K £ (1 - cos0S) - ko (1 - cos0 o ) = \ (((A* - (£1 - £2)) c) 2 - (A*c) 2 ) 



Therefore, 



so 



(12) 

These calculations lead to the first important subsidiary distribution for quantifying the OSMMTS 



error. 



Definition 9 Let n = n (tt^tS, £1,62,63) = cos 05 - cos0 o . Ffcen 7/ is said £o have an M 0 (tl,t$,tl) 
distribution. 



Hence, 



and 



1 - cos 0 O = 1 - cos 0S + V 



' T/f c^4«S + (/f c) 2 



= < 



± (/1 - (ei - ca)) c 



2 ( rf?-((/r-(^i-^))c) 2 

y (l-COS^+T?) 

+((/r-(^i-^))c) : 



If ^ = (/* - ( £l - e 2 )) c ~ iV (tfc, 2c 2 <t 2 ) , then we have 



(41 ~ sll) - (4 - sl x ) = T/fc 



(1 - cos 0£) 

Definition 10 Wnen £ ~ AT (tfc, 2c 2 a 2 ) and n ~ M 0 (fl,*5, tg), tfien 



1 



d? - (A*c) : 



+ (Ac) 2 ±6j2 



d?-£ 2 



cos 05 + 77) 



+4? 



(13) 



=F/rc 



2(y 2 -2/i) I 1 ^ \(l-cos^) 



2l 4^Li + (A* c) 2 ±^ 



/T7 

'Y V(l-cos0S 



+ »?) 



ts distributed according to the Mi (t^t^.ts) probability distribution function. 

We also have that when & = (/ 2 * - On ~£3))c " N (/ 2 *c,2c 2 a 2 ) and t; " M 0 (tl.t2.t3). then 



2 (y3 - yi) 



N 



2 I j ~ {n t \ + (/ 2 'c) 2 ± 6 



til—"! 

COS Oq + T)) J 



+ £ " Mi(*l,*2,*8) 



(l-cos05)y y v^-^^o 

Definition 11 //X " Mi (*;,*£,*£) and y " Afi (tJ.tJ.tS), and X and 7 are independent, then 

X-Y~ M 2 (tlt*,tl) 

Definition 12 // J " M 2 (tj , tg) and K ~ M 2 (t;,t5,*S), and J and are independent, then 

(mi _ m2 ) J 2 + ^1 - ^ tf 2 " M 3 (tl.ij.tj) 
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By (13), we have 

2(y 2 - yi ) M ^M) 

with a similar statement for Ifii^hjgZfiiJ - m x (t^t^q). Let 

\mi — 7712 \mi — 77Z2 / / 

with similar results for other forms of (x3>2/o) corresponding to (xo>yo) ? where 

y = fSj -*ii + (*l-ai) + (y2-y?) 

1 2(2/2-2/0 

2 2(y 3 -yi) 

Note how the value of mi and 7712 are unaffected by consideration of the t*. 
Then 

b * = g2i-gii + (g2-gi) + (y2-yi) (S21-gll)-(*21-*ll) 

1 2(2/2-2/1) 2(2/2-2/1) 

,* _ S 22 ~ 5 12 + ( x 3 ~ x l) + (yl ~ V\ ) . { s 22 ~ ~ ( 5 22 ~ s 12) 

°2 ~ ~ v r 



or 



2(y3-yi) 2(2/3 -yi) 

( 5 21 ~ S ll) - ( S 21 ~ S ll) 



6* = 61 + 



62 = 62 + 



By Definition 10, this means 



2(2/2 -yi) 

{ S 22 ~ S li) ~ ( 5 22 ~ 5 12) 

2(ya - yi) 



Therefore, 



and by Definition 11, we have 
Likewise, 



61 -61 - A#i(t;,«s,t5) 
63-62 " Mi («;,*;, ^) 

a* _ a; 0 = &2=6l_ 

u " mi— 7712 till— T7l2 

_ oi-tn-(ta-ti) 

(mi - m 2 ) (xq - x 0 ) ' Mi {t\,t* 2 , t$) 

w-» = m '(S) + ^-( m >(S) +6 ') 

m . (J£^L _ 6,-6, \ 

A ymj— 1712 TTll — T7l2 / 

=fe «*S " 62) - (61 - 60) 



and by Definition 11, we have 



30. The Error Likelihood Ellipse 

Definition 13 The Error Likelihood Ellipse (ELE) for the OSMMTS position (x 0 ,y 0 ) is 

( ig -x 0 \ 2 ( yZ-yp \ 2 _ 

[tz%) + {w^r 2 ) - R 

with (mi - m 2 ) (arj - x 0 ) " M 2 (t^t^t^) and ^1 - (2/^ - yo) " M 2 am * w/iere tfie constant 

R depends on the two M 2 J ,^2^3) distributions. 



t 
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However, we have 

p — ( z n- x o i i ( vq-vo \ 2 

m l 

(m 1 -m 2 )(x;-x 0 ) 2 +(l-^)(y;-y 0 ) 2 

(mi-m 2 )(l-^) 
( mi - m2 )( g ;^xo) a +(l-^)(y;-y 0 ) 3 

/ (m 1 -m 2 ) 2 A 



or 



K - m 2 ) W - x 0 ) 2 + (l - ^) (Vo - J/o) 2 = {mi ^ 2)2 R 
By Definitions 11 and 12, we have 

77li 

Hence, the likelihood of the Target T being within the Error Likelihood Ellipse with constant < mi ~™ 2)2 R 
is calculated by the probabilities of M 3 (£{,^,£3). 



I 



PQIC Documentation — OSMMTS Patent Application Technical Documentation 33 

31. Error Likelihood Ellipse Calculation 

The Error Likelihood Ellipse is the totality of Target Positions Reports for a particular set of three SDU's 
produced by varying the observed arrival time data according to an assigned error distribution. The resulting 
* (x, y) position that has maximum likelihood within this ellipse is judged to be the most likely "real" position 
of the target. The same algorithms used to calculate an individual Target Position Report shall be used to 
calculate the Error Likelihood Ellipse. 

32. Error Likelihood Ellipse Standard Methodology 

The PQIC OSMMTS Error Likelihood Ellipse Standard Methodology is to assume the arrival time data 
varies according to a normal distribution of mean zero and standard deviation a, where a is either estimated 
from calibration data or the result of a Bayesian estimation process. This means the observed arrival time 
data is always accepted as the center of the error distribution. 

The algorithms found below implement such an error distribution. However, such algorithms may be 
amended at any time to more accurately reflect the true error distribution in any particular PQIC OSMMTS 
System implementation. 

33. Error Likelihood Ellipse MAPLE Algorithm 

Given the position information of the SDU's that receive the arrival time data, the following MAPLE 
algorithm calculates the Error Likelihood Ellipse based on the normal error distribution affected arrival 
times. These algorithms call the Target Position Report algorithms by reference, so that the Error Likelihood 
Ellipse algorithms are independent of the particular methods found in the Target Position Report code. 

The M 0 distribution and the Target Position Report algorithms are needed to calculate the Error 
Likelihood Ellipse. All other M distributions are only needed when quantifying the bounds of the error 
ellipse, or for making multivariate calculations on multiple Target Position Reports and Error Likelihood 
Ellipses. 

34. ( ELE Calculation 34 ) = 
MOGEN: = proc(tl ,12 ,t3 ,a,6,r) 
local ans, el , e2, e3,jl i j2 1 j3 1 idx 1 tx\ 

options' Copyright2Q03?QI Consulting All Rights Reserved 1 ; 
de5cr^p^i^?n ,, 0SMMTSuM0uGeneration u Calculation ,, ; 
idx: = 0; 
with (stats)] 

for el fromatob do jl : = statevalf [cdf , normald [0,0.1]] (el /r) — $tatevalf[cdf, normald[0 } 0A]]((el -l)/r)] 
for e2 fromatob do j2: = statevalf [cdf , normald[0 i 0A])(e2 /r)-statevalf[cdf y normald [0, 0. 1]]((e2 -l)/r); 
for e3 fromatob do j3 ': = statevalf [cdf , norma Id [0,0.1] ](e3/r) - statevalf [cdf , normald [0 , 0. 1]]((e3 -I)/ r)] 
tx: = evalf (osmmts(tl + el /r,t2 + e2/r,t3 + e3/r)[],5); 
if (type(tx[3\, undefined)) then 
else 

if (abs(tx[A)) > 1) then 
else 

idx: = idx + 1; 

ans[idx]: = tx[2]Jl * j2 *j3)\ 

end if ; 

end if ; 

end do ; 

end do ; 

end do ; 

return ans y idx\ 

endproc) 



t 



34 PQIC Documentation — OSMMTS Patent Application Technical Documentation 

J: = proc(dt,cnt,msh) 

local idx , inx , iny , new , smm , minx , maxx , miny , maxy ,g,gx,gy,M\ 
options ' Copyright 2003PQ I Consulting All Rights Reserved 1 ; 
description n OSMMTS u MO u Processing" ; 
minx: = 0; maxx: = 0; 
miny: = 0; maxy: — 0; 
=0; 

for idx fromltocnt do new: = d*[idx]; 
inx: = floor(msh * new [1]); 
iny: = floor (msh * new [2]); 
if (a55^ned(M[ma: ) iny])) iften 
M[inx, iny]: = M[mrr, iny] + new [3]; 
else M[inx,iny]: = new[3]; 
end if ; 

if (inx > maxx) then 
maxx: = inx; 
end if ; 

if (inx < minx) then 
minx: = inx; 
end if ; 

if (iny > maxy) then 
maxy: = inx; 
end if ; 

if (iny < miny) then 
miny: = inx; 
end if ; 
end do ; 

for idx from minx to maxx do for iny from miny to maxy do 

if (assigned (M[idx, iny])) then 

if (M[idx,iny] > g) then 

g: = M[idx, iny]; 

#x: = idx; 

#y: = iny; 

end if ; 

else M[idx,iny]: = 0; 
end if ; 
end do ; 
end do ; 

return evaZ (M) , minx , maxx , miny , maxy ,gx,gy,g\ 
endproc; 

This code is used in chunk 12. 

35. Error Likelihood Ellipse Example #1 

36. ( ELE Example #1 SDU Coordinates 36 ) = 
(xUyl): = (2,3); (x2,y2): = (4,2); (x3,y3): = (3,1); 

This code is used in chunk 12. 

37. Error Likelihood Ellipse MAPLE Calculation For Example #1 

> with(stats): 
^ mesh * — 3 5 * 

> TPR : = (evalf (osmmts ( 16/10 , 2 , 175/100) [1 . . 2] , 10) ) ; 
TPR := 2.547503694, 2.097308261 
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> RES : =M0GEN (16/10,2, 175/ 100 , - 12 , 12 , mesh) : 

Warning, these names have been redefined: anova, describe, fit, 
importdata, random, statevalf, statplots, transform 

> MZ:=J(RES[1] ,RES[2] ,mesh) : 

38. Error Likelihood Ellipse Example #2 

39. ( ELE Example #2 SDU Coordinates 39 > = 
(xl,yl): = (3,2); (x2,y2): = (2,4); (x3 y y3): = (1,3); 

This code is used in chunk 12. 

40. Error Likelihood Ellipse MAPLE Calculation For Example #2 

> with(stats) : 

> mesh: =35; 

> TPR : = (evalf (osmmts ( 16/10 , 2 , 175/100) [1 . . 2] , 10) ) ; 
TPR := 2.097308261, 2.547503694 

> RES : =M0GEN ( 16/ 10 , 2 , 175/ 100 , - 12 , 12 ,mesh) : 

Warning, these names have been redefined: anova, describe, fit, 
importdata, random, statevalf, statplots, transform 

> MZ:=J(RES[1] ,RES[2] ,mesh) : 
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41. Graphical Methods For Target Position Report And Error Likelihood Ellipse Display 

The analytical methods demonstrated here relate the graphical depiction of the Error Likelihood Ellipse to 
the relative position of a given Target Position Report and the most likely position value within the Error 
Likelihood Ellipse. These depictions are not meant to list an exhaustive summary of all such graphical 
methods. The policy of the PQIC OSMMTS System is to leave such methods to the particular needs of an 
implementation's circumstances. The methods found herein are for demonstration purposes only. 

> with (plots) : 

Warning, the name changecoords has been redefined 

> F : =PL0T3D (POLYGONS (seq( seq( [ 

> [(i-l/2)/mesh, (j-l/2)/mesh,MZ[l] [i,j]] > 

> C(i+l/2)/mesh, (j-l/2)/mesh,MZ [1] [i,j]] > 

> [(i+1/2) /mesh, (j+1/2) /mesh, MZ[1] [i,j]] , 

> [(i-l/2)/mesh,(j+l/2)/mesh,MZ[l][i,j]]] , 

> i=MZ[2] . .MZ[3]),j=MZ[4] . .MZ[5])), 

> COLOR(ZHUE) , AXESSTYLE (NORMAL) ,AXESLABELS(' ! X n , "Y" , "Z") , 

> VIEW (TPR [l]-l/2. .TPR[l]+l/2,TPR[2]-l/2. . TPR [2] +1/2, DEFAULT) ) : 

> SI : =PL0T3D (POINTS ( [MZ [6 . . 7] /mesh , 0] ) , COLOR (RGB ,0.0,0.0,1.0), SYMBOL (CR0 

> SS,50)): 

> S2 : =PL0T3D (PO INTS ( [TPR [ 1 . . 2] , 0] ) , COLOR (RGB ,1.0,0.0,0.0), SYMBOL (CROSS , 5 

> 0)): 

> display({F,Sl,S2}); 




> F: =PL0T3D (POLYGONS (seq(seq([ 

> [(i-l/2)/mesh,(j-l/2)/mesh,MZ[l] [i-l,j-l]] , 

> [(i+1/2) /mesh, (j-1/2) /mesh, MZ[1] [i ,j-l]] , 

> [(i+1/2) /mesh, (j+1/2) /mesh, MZ[1] [i, j]] , 

> [(i-l/2)/mesh,(j+l/2)/mesh,MZ[l] [i-l,j]]] , 

> i=MZ[2]+l. .MZ[3]),j=MZ[4]+l. .MZ[5])), 

> COLOR(ZHUE) , AXESSTYLE (NORMAL) , AXESL ABELS ("X" , "Y" , "Z") , 

> VIEW (TPR [1] - 1/2 .. TPR [1] +1/2 , TPR [2] -1/2 . . TPR [2] +1/2 , DEFAULT) ) : 

> S : =PL0T3D (P0 INTS ( [TPR [ 1] , TPR [2] , TPR [3] ] ) , COLOR (RGB ,1.0,0.0,0.0), SYMBOL 

> (BOX, 50)): 

> display({F,S}); 
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42. The Target Position Report Vector 

The Target Position Report Vector (TPR) is the quadruple defined by the Calculated Target Position 
(rr 0 ,2/o), the Maximum Likelihood Target Position {x^y^), associated error likelihood ellipse ^ f^sp j + 

{£=^) 2 = R > ^ d the associated likelihood value for ( m *-^) a fl ~ M 3 (ti.tS.tJ): 

TPRV= ^ 0)J/0 ), K) ^)^^y + ^A^y = R ,^l^t R - Jf, j 

The value of R is chosen so that ( mi ~™ 2 ^ R is at the 100 (1 - a) % percentile of the M 3 (*i,<5,*|) 
distribution. 
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43. Containment Policies 

44. Definitions and Remarks 

Definition 14 A TPR is said to be accurate if the actual position of the target is inside the associated 
ELE. The value of R determines the likelihood of this event. 

Definition 15 Any calculation algorithm used to produce a set of numerical values intermediate and inferior 
to the TPR is called an analytical step. 

Definition 16 An analytical step is called a mitigation if it is taken before the time data {ti, £2, .} 
are collected. 

Definition 17 An analytical step is called a optimization if it occurs after the time data {t\, t<i, .. . , tk> . . .} 
are collected. 

Remark 18 The purpose of mitigation steps is to reduce the error variance a. 

Remark 19 The purpose of optimization steps is to increase the likelihood of an accurate TPR. 

Definition 20 An irregularly occurring y non- analytical step taken at any time to accomplish the same goals 
as mitigation and optimization is called ad-hoc. 

Definition 21 The collection of ad-hoc, mitigation, or optimization steps taken in an implementation 0} the 
OSMMTS is called the system's containment policies, and referred to individually as a system containment 
policy. 

45. The Demerit System 

The OSMMTS Demerit System is an ad-hoc containment policy that acts simultaneously as a mitigation 
and an optimization. Under this system, the three SDU's chosen to calculate the TPR are those three that 
are most likely to produce the "best" TPR based on past performance (thereby making it an optimization 
step), by way of reducing the variability of the utilized data (thereby making it a mitigation step). 

Suppose there are n-many SDU's, however, only k < n many receive a signal within the reception 
window. There are (£)-many combinations of SDU's, and (*)-many combinations of the A>many that 
receive the signal taken three at a time. Each SDU has three values associated with it at the beginning of 
each processing cycle, namely its non-negative Demerit Count, its positive History Total, and its possibly 
null Boolean Confirmation Value. At the beginning of all processing, the demerit count for each SDU will 
be zero, the history total with be one, and the confirmation value will be null. The confirmation value at 
the beginning of the processing cycle is determined by its observed value during the confirmation cycle. At 
the end of a processing cycle, the demerit count and history total are determined by the steps below, and 
the confirmation value is set back to null. 

For each processing cycle, and for each of the (^-many combinations, the following steps determine the 
end-of-processing-cycle demerit counts and history totals. 

1. Set the likelihood value A. 

2. Eliminate those ri-many combinations that are collinear. 

3. Eliminate those T2-many combinations that do not all have positive history totals and TRUE confirma- 
tion values. 

The SDU's involved in the (r\ -f r 2 )-many combinations eliminated in Steps 2-3 are called deficient 
for the current processing cycle. This designation is removed at the beginning of a new processing cycle. 

4. Among the remaining, i.e., qualifying combinations, choose the combination of three SDU that collec- 
tively have the minimal sum of demerits. 

5. In case of a tie in Step 4, use the combination with the largest history sum. In case of a further tie, 
choose the combination with the smallest individual demerit count. In case of a last tie, randomly 
choose uniformly among the finalists. 

The combination so chosen is called the calculating combination, and the SDU's involved are called 
the elected SDU's. Increment the history total by 1 for each elected SDU. 

6. Subtract two demerits from the count for each elected SDU. Recall the demerit count for an SDU cannot 
become negative. 
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7. Calculate the TPR using the calculating combination. 

8. Calculate the A-ELE for the calculating combination. 

9. Calculate the ^(3) - (ti + T2)^-many TPR for all other qualifying combinations. Each of these TPR 
is called an Alternate Position Report (APR). 

10. For each APR calculated in Step 7, if the APR falls outside the A-ELE, then add one demerit to the 
count for each SDU involved in the APR. 

11. For each APR calculated in Step 7, if the APR falls inside or on the A-ELE, then subtract one demerit 
to the count for each SDU involved in the APR. Recall the demerit count for an SDU cannot become 
negative. 

12. Add one demerit for each SDU that does not report a positive confirmation. 

13. When the demerit count for an SDU exceeds the Warning Threshold, send an alert to report a frequently 
deficient SDU. 

14. When the demerit count for an SDU exceeds the Terminal Threshold, shut down communication with 
the SDU and do not consider it further (by setting its history total to zero) until explicitly reset. Also 
send an alert to report a failed SDU. 

15. These steps are in addition to the disabling of an SDU if proper query responses are not confirmed 
during the receive phase. 
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46. Mitigations For Reflections 

Given aTPR {x 0i y 0 ) and multiple detection times at an SDU, say {£1,^2, £3, . . .}, where t\ is used to generate 
(xq, t/o) , the question of mitigating the error in the TPR due to reflections 13 of the "true" signal may be 
stated as follows: What is the earliest time a reflection of the "true" signal could be detected at an SDU, 
and what is the latest such time. Any signal detected between these two extremes would be considered a 
reflection, and therefore not used in the calculation of a TPR. Any signal detected "too early" or "too late" 
would indicate the detection time data is corrupt and should not be used to calculate the TPR. 

Axiom 22 // the signal detection times recorded during the Receive Phase are inconsistent with a calculated 
TPR, the TPR should not be used. 



47. Definitions 




Figure 10: Generic Straight Line Obstruction 



Consider a generic straight line obstruction as depicted in Figure 10. Even though the following develop- 
ment requires a straight line obstruction, a more generalized parameterized formulation may be substituted 
for 14, such as a Bezier curve or trigonometric function. 

As s goes from 0 to 1, point (x 8} y 8 ) goes from L = to R = (x r ,2/ r ), so that 

(x 8i y 8 ) = s(x r ,y r )'r{l-s)(x h y l ) ^ 
= (xi + s (x r - xi) , yi + s (y r - yi)) 

Then given SDU position (xi,yi), TPR position (xo>2/o), and obstruction endpoints L and R t we wish 
to find 



and 



subject to 



rfmin = nrin {d\ + d 2 } 

0<«<1 



4ax = max {di + d 2 } 

0<8<l 



d 2 = d\ + d\ - 2d x d 2 cos (it - 2a) 



13 



The term "reflections" also includes signal behavior characterized by "echos," "secondary," and "phantom." 
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since the angle of incidence a would be consistent "inbound" to, and "outbound" from, the obstruction. 
Define the Window for Reflections Between I and T to be 

dmm < t < dmax (15) 

c c 

Note that at point (x Si y s ) we have 

{x 0 - x s ) * (x r - x 8 ) + (y 0 ~ 2/5) * (2/r - Vs) 



a = arccos 



y (2* - z fi ) 2 + (y 0 - y 8 ) 2 \/{x r - x s ) 2 + (y r - y s ) 



so that 

cos (7r — 2a) = — cos 2a 

= 1-2 cos 2 a 

_ _ 2((xo-x 8 )*(x r -x s )+(yo-y a )*(y r -y a )) 2 
((xo-x a ) 2 +{yo-y a ) 2 )((x r ~x 3 ) 2 +(yr-ys) 2 ) 

((X 0 ~Xi)-S (X r - Xi)) * ((1 - 5) (x r ~ Xl)) 

1 o v + ((go -yi)- s (vr -m)) * ((1 - g) far - y/)) , 

/ ((X 0 - X,) - S (x r - X/)) 2 \ / ((1 - s) (X r ~ Xl)) 2 \ 

\ + (fao -Vi)-s (y r - y/)) 2 / V + ((1 - s) (y r - y t )) 2 J 

= 7 00 

and 

di (s) = y ((an - x t ) - s (x r - xi)) + ((yi - y/) - s (y r - y/)) 

and 



d 2 (s) = \/ ({x Q -x t )-s {x r - xi)) 2 + (fao ~yi)-s {y r - yi))' 
which means d 2 (5) may be viewed as a function of s. 



48. Policies 

The Mitigation for Reflections Method calculates 

^ = \ min {di + d 2 } and ^ = \ min {d x + d 2 } 

subject to 

d 2 (s) = d? (5) +<%(s)- 2dx (s) d 2 (5) 7 (s) 

for each combination of SDU position, possible target position, and L and R positions for a straight line 
obstruction that would potentially corrupt the TPR through the detection of reflections (see again Figure 
10). 

Given a TPR (x 0 , yo) and multiple detection times {t\ , t 2 ,t 3 , . . .} at SDU I, where ti is used to calculate 
the TPR, if 

dm in ^ . - d m ax 

— S *i s — 

c c 

for a// i = 2, 3, . . ., then the TPR is accepted, i.e., there is no need for a mitigation for reflections 
If, however, 

d • d 

U < or < U 

c c 

for any i = 2, 3, . . ., then the TPR is rejected, and a new TPR is calculated using the next highest priority 
qualifying combination of SDU's. This cycle continues until a TPR is accepted or the qualifying combination 
of SDU's is depleted. 



49. MAPLE Implementation Code 

50. (Mitigations For Reflections 50 ) = 
( Mitigation Basis Algorithm 51 ) 

( Mitigation Extrema Algorithm 52 ) 

This code is used in chunk 12. 
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51. (Mitigation Basis Algorithm 51 ) = 
mitig: = proc(xx.O , yy.O) 

locals, x_s , y.s , d, dA , d.2 ,g,sO, xO.s , yO.s ; 
global x-l , yJ , x.r , y.r , x_i , y.l ; 

options ' Copyright 2003PQI Consulting All Rights Reserved' ; 

description " OSMMTS u Mit igat ion u Algor ithm" ; 

if (xx-0 , yy_0) = (x_i , y-1 )£/ienreturn -1, -1, — 1, -lend if ; 

x.s: = xJ + 5 * (x_r — xJ); 

y_s: = yJ + s*(y.r-yJ); 

dA : = sgrt ((((x_i - xJ) -s* (x.r - x.l)) © 2) + (((y_i - yJ) - s * (y_r - y.l)) © 2)); 
d_2: = sqrt((((xx.O - xJ) - 5 * (x.r - xJ)) © 2) + (((yy.0 - yd) -s* (y.r - y_/)) © 2)); 
d: - sqrt((xA - xx_0) © 2 + (y_l - yy.O) © 2); 

5: = simplify ( 1 - 2 * (((((xx_0 — xJ) — s * (x_r — x_/)) * ((1 - 5) * (x.r — xJ)) 
+((yyj?-yj)-s*(y_r-yj))*((l-s)*^ 

-5 * (y.r - yJ)) © 2) * (((1 - s) * (x.r - xJ)) © 2 + ((1 - s) * (y.r - yJ)) © 2)) ) ) ; 
sO: = fsolve(dA © 2 + d_2 © 2 - 2 * dA * d.2 * £ = d,s,s = 0..1); 
if substring (convert (sO , string), I.. 3) <> "f so" then 
return sO,subs(s = sO,x.s),subs(s = sO,ys),subs(s = sO,dA + d-2); 
else return —1, minimize (#, 5 = 0..1), maximize(g,s = 0..1), —1; 
end if ; 

This code is used in chunk 50. 

52. (Mitigation Extrema Algorithm 52) = 
xmit: =proc() 

local res , minres , maxres , u, v; 
global x_0, y_0; 

options 1 Copyright 2003PQI Consulting All Rights Reserved 1 ; 
description n OSMMTS u Mitigation u Extrema u Calculat ion"; 
minres: = mitig (x.O , y_0)[4]; 
maxres: = minres; 

for u Trom — 30£o30 do for v from — 30£o30 do res: = mitig(x.O + u/10, y_0 + v/10)[4]; 

if res <> —lthen if res < minres then minres: = res end if ; 

if res > maxres then maxres: = res end if ; 

end if ; 

end do ; 

end do ; 

return minres, maxres; 
endproc; 
This code is used in chunk 50. 

53. Sample Calculations And Results 

> (x.l,y.l):=(l,5); (x.r ,y_r) :=(2,2) ; (x_0,y_0) :=(4,3) ; 

> (x_l,y_l):=(3,6); 

xJ, yJ :— ' 1, 5 
x.r, y.r := 2, 2 
x.0, y.O := 4, 3 
xA, yA := 3, 6 

> mitig(x_0-l ,y_0+3) ; 

-1,-1,-1,-1 

> xmit(); 

3.493116110, 9.112507778 
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54. Optimal Placement Of Detection Units 

A 



! V ! 

! • i 

77/ i'' N i 7K 



Figure 11: Poor Placement Of SDU Units 

The optimal placement of SDU's depends on the OSMMTS goals for the implementation under consid- 
eration. In general, a maximum likelihood approach to SDU placement works well for many applications. 
Its goals address the probability of an accurate TPR calculation with the first qualifying set of SDU's within 
a particular SDU arrangement. The probability of such placement is maximized among all convex polygons 
that encloses all possible target locations, or those target locations that are of interest in the implementation, 
or where the addition of any more SDU's would not enclose a significantly larger area than is covered by 
smaller set of SDU's. Such an approach uses the following rules for calculating a set of exact SDU positions 
for optimal use. 

Optimized SDU Placement Rules 

1. SDU's should be placed where the connected units form a convex hull (see Figure 12 for an example of 
a convex hull, and Figure 11 for an example of a non-convex hull). 

2. The convex hull should cover as much of the possible target locations as possible. 

3. SDU's are placed so that no three units are collinear. 

4. SDU's are placed so that no two units have the same x-value or y-value, as viewed in a coordinate grid. 

5. SDU's should be placed to maximize the likelihood that d{ > \fc\ c, for all i. 

No consideration is made in these rules for obstructions, or other exceptions in a clear area of target 
detection. Under the OSMMTS, such considerations are handled by mitigations. Furthermore, if the set of 
possible target locations, or those target locations that are of interest in the implementation, contain "gaps," 
regardless of size and shape, such gaps shall not be considered in the optimal placement of detection units, 
since their consideration may be completely addressed by the methods found in the containment policies. 
Specialized implementation of the OSMMTS may introduce special SDU placement considerations, which 
will be documented under separate cover in the particular implementation documentation. 

Figure 13 shows an arrangement of five SDU's that satisfy the guidelines embedded in the placement 
rules. Since the position of the origin is arbitrary, the coordinate center of the region may be translated to 
any local coordinate system within any particular OSMMTS implementation. 
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Figure 12: Better Placement Of SDU Units 




Figure 13: Best Placement Of SDU Units 
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55. The Convex Hull Principle 

Since the area of any segmented convex hull is maximized when the distances between the vertices of the 
segments are equal 14 , the "best" convex hull, i.e., that which follows the placement rules, is based on a circle 

14 Since all convex polyhedra may be subdivided into triangles, it suffices to prove the following claim to make this assertion. 
Claim 23 Given a fixed perimeter length, the triangle with maximum area is equilateral. 

Proof. Let a > 0, b > 0, and c > 0 be the length of the sides of a triangle. Then the area of such a triangle is given by 

A = y/s (s - a) (s — b) (s — c) 

where s = | (a + b + c). It remains to show that A achieves its maximum when o = 6 = c, given a fixed value of a 4- 6 + c. 
Consider the Lagrange Multiplier Function K: 

K(a,b,c) = A- \(a + b + c) 



Then 



where 



Bi = 



Similarly, we have 



where 



Setting 

we have 
or 

or 

which means 
Furthermore, 

gives 

or 

which means 
In fact, 



dK dA _ Bi 
da ~ da ~ AA 

(-a + 6 + c) (a — 6 + c) (a + 6 - c) 
- (a + b + c) (a - 6 + c) (a + 6 - c) 
+ (a + b + c) (-a + b + c) (a + b - c) 
+ (a + 6 + c) (-a + 6 + c) (a - 6 + c) 

— - — - A- — - A 
db ~ db ~ AA 

(-a + 6 + c) (a - 6 + c) (a + 6 - c) 
, +(a + 6 + c)(a-6+c)(a + &-c) 
** 2- l -(a + 6 + c)(-a + 6 + c)(a + 6-c) 
+ (a + 6 + c) (-a + 6 + c) (a - 6 + c) 



dK _ dK 
da ~ db 

AA AA 



f _ ( a + + c ) ( a - b + c) (a + 6 - c) If + (a + 6 + c) (a - 6 + c) (a + 6 - c) 1 
| + ( a + 6 + c )(_ a + 6 + c )(a + fe-c) J \ - (a + 6 + c) (-a + 6 + c) (a + 6 - c) J 



a-6-j-c=-o + 6 + c 
a = 6 



dK _ Q dK 
da ~ ~ dc 



/ -(a + 6 + c)(a-6 + c)(a + 6-c) \ _ f + (a + 6 + c) (a - 6 + c) (a + 6 - c) 1 
\ + (a + 6-i-c)(-a + & + c)(a-& + c) J-\ -(o + 6 + c)(-a + b + c)(a-6 + c) J 



a+6-c= — a + 6 + c 



a = c 



a 2 K 4ABu - ^ 
0a 2 ~ 16yl 2 



where 



Bn = < 

Then a = b — c gives 



-2 (a - 
+2( 



: (a - 6 + c) (a + 6 - c) "| 

+4a(-a + 6 + c) > 

a + 6+c)(-3a + 6 + c) J 

- { 5 } = 



and 



«— - (I*-) (S~) - J ^ - ^ 
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(which would be an "infinite-many sided polygon"), and the n-many SDU's would be at positions 

re- = r(^— j +t sm(— Jj 

for k — 0, 1,2, ... ,n — 1, and where r is determined by the effective range of the possible target positions, 
relative to an arbitrarily chosen "origin" for location purposes only. Note how such an arrangement au- 
tomatically satisfies Placement Rules #1 and #3, and Placement Rule #5 is violated only when a target 
comes between the SDU's on the "edges" of the containing circle. This consideration may be addressed 
through the choice for r, which also covers Placement Rule #2. 

The value of r is called the Radius of Coverage, and the placement of the units so described is called 
the SDU Arrangement. 

For example, for n = 8, the SDU's would be at positions eT 1 , for k = 0, 1, 2, . . . , 7. This translates to 

r (cos n + i sin n) 

where 

7T 7T 37T 57T 37T 77T 

*~ ' 4' 2'T ,7r 'T'T'T 

(see Figure 13). 

Note how, for n = 8, there is a problem with having SDU's at both 



/ 7T . 7T\ / 3?T 3?T 

^cos — + 1 sin —j and r I cos — + 1 sin — 



since sin J = sin^, which violates Placement Rule #4. However, choosing any prime n > 4 satisfies 15 
Placement Rule #4. 

The optimal placement of SDU's may now be characterized analytically in terms of the "point of 
diminishing returns." 

Condition 24 Given a coverage function A(n,r) that describes the area covered by the n-many SDU ar- 
rangement with radius of coverage r } find 

n* (r | 0 = m|n (n,r) < £ j 

Then given r, n* (r) is the minimum number of SDU's necessary to satisfy the placement rules, and 
where the coverage area instantaneous rate of increase first drops below f units of area. By choosing a 
sufficiently robust value for f , the optimal number n m may be explicitly calculated. 



56. Formulation of A (n, r) 

Since the SDU arrangement is completely symmetric for any two SDU's and the center of the containing 
circle, consider the sector spanned by two SDU's A and B } and the center of the containing circle O (see 
Figure 14). Note that the position and orientation of this figure is arbitrary, and the results of this section 
do not depend on the particulars of such aspects. 

Lemma 25 The area of AAOB is \r 2 sin (it (l - 



and 

3\2 



BJK ^ ( - 4 " 2) ' ffi (-<^-a§)° 4 . 

a=6=c 16 (^I a 2) 2 3a* 3 V VZJ 

regardless of the value of a. Similar calculations confirm < 0 and | < 0 regardless of the value of b or 



c. 



15 If n is prime, then 

W\ 5* W2 => cos 

and 

wi W2 => sin 

since n would not be divisible by 2, or by any W{. Hence, no two SDU's would have the same x- value or same y- value. 



a=b= 


c 




2wi fi- 


^ COS 






ll 


n 


2w\7T 


sin 




n 


n 
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B 



Figure 14: Optimal Placement Sector 



Proof. In Figure 15, ACOB is a right triangle with height OC and base CB. Since OB = OA, then 
/? = corresponding to ZOAC and ZOBC. 
Then 



which means 



OC = r sin/? and CB = r cos/3 



Area Of ACO£ = -r 2 sin 0 cos 0 

it 



Since ACOB is similar to ACOA, then 



Area Of AAOB = r 2 sin /J cos £ = -r 2 sin 2/3 

_ 2tt 



In the SDU arrangement, we have a = for n > 4. Hence, the area covered by the convex hull of 

l n - - 

2 

-2 



SDU's is incrementally given by \r 2 sin 2/3, where 



P 2 2n 



=7r (V) =7r G-^) 

Hence, for an n-many SDU arrangement, we have 

Area Of AAOB = ^r 2 sin 2/3 = ^r 2 sin ^tt ^1 - J 

■ 

Claim 26 A (n,r) = § r 2 sin (ir (l - £)) 

Proo/. The coverage region is comprised on n-many exactly sized sectors of the type found in Figure 
14. The result follows immediately by Lemma 25. ■ 

For a fixed r, as n -> oo, we have 

and 

sin K 1 -^))- >0 
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In fact, 



c 



B 



/ 

/ 

/ 

,'r 



o 



Figure 15: Sector Division 



n— >oo 



.2 sin(^(l-^)) 



lim §r 2 sin(7r(l-£)) = r 2 Urn 

— 2 v v nn n-yoo 

= r 2 lim V J_ /; 

n— ►co ~ T7 

= -7rr 2 lim cos (V (l - £)) 



n— kx> 



= -7rr 2 (-l) 

.2 



= 7rr 



which is the area of a perfect circle of radius r. 
Now by Claim 26, we have 



dA 
dn 



(n,r) = r 2 (±sin (l - ^)) + ^cos (l - 2))) 



the graph of which is contained in Figure 16. This may be used to calculate n* for a given £. 
Finally, solving A(n,r) for r gives 



= \l 2 n ACSC H 1 - 2 n)) 



and if the approximate coverage area A and the approximate radius of coverage r axe given, then define 



n„ = esc /zj + 
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Figure 16: in-Many Vertex Convex Area Of Radius r 



to be the optimal number of SDU's, where z — \x satisfies 16 

A (z — 7r) + 7rr 2 sin z = 0 



57. MAPLE Implementation Code 
> X : = (n/2) *r ~2*sin(Pi* ( 1- (2/n) ) ) ; 



> solve(A=X,r) [1] ; 



1 2 

X:=-nr 2 sin(7r(l--)) 
I n 



V2,/nsin(^^M 
V n 

-.:_/ *("- 2 ) l 

n sin( j 



n 



> subs ( A=1000 , n=9 , A=X) ; 



> f solve C/,,r,0. . 100) ; 



1000 = ^r 2 sin(^) 



18.59345062 



16 



Note that if A {n — 7r) + irr 2 sin /x = 0, then 

2 



fr 2 sin(7r(l-£))| 2a 
2 v v nil in=-4j-yi csc 



^ ^ CSC _ 2 



= i 4csc/isin(7r(l-^ r )) 

= A esc /x sin ((7r — \itr 2 sin/x)) 

= /Icsc/isin ((7r + (^i — 7r))) 

= A esc /x sin \i 

= ,4 
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> solve ( A=X, n) ; 



2A 



> 
> 
> 



r 2 sin(RootOf {A -Z + irr 2 sin(_Z) -Air)) 
subs(A=1000,r=18.25,A=X) ; 

1000 = 166.5312500 nsin(7r(l - -)) 

n 

f solve('/.,n,4. .100); 

12.10617840 

plot 3d (X , r=0 . . 10 , n=3 . . 20 , 

or ientat ion= [55 , 60] , axes=normal , grid= [10 , 20] 

) ; 




> dXdn:=diff (X,n); 



i 9 r 2 cos(7r (1 

dXdn := - r 2 sin(7r (1 - -)) + 

Inn 



n 



> 
> 
> 



plot3d(dXdn,r=0. . 100, n=3. .50, 

or ientat ion= [25 , 60] ,axes=normal ,grid= [10 

,20]); 



> C : =subs (r =25 , dXdn) ; 



625 2 625cos(tt(1 

C:=^sin(7r(l-£)) + 

In n 



n 



> solve (C=10,n) ; 



-2 7r/(RootOf((15625sin(_Z) 2 - 15625) -Z 2 + (-31250 tt sin(_Z) 2 + 31250 tt) ,Z 
+ 15625 7r 2 sin(_Z) 2 - 15625 tt 2 + 15625 sin(.Z) 2 - 1000sin(.Z) + 16) - tt) 
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58. (Coverage Area Optimization 58) = /* Embedded In MAPLE Code*/ 

This code is used in chunk 12. 



59. Sample Calculations And Results 

> f solve (C=10,n=3. .50) ; 

13.62494565 

> f solve (C=5,n=3. .50) ; 

17.21207895 

> f solve (C=l,n=3. .50); 

29.51887628 
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( Mitigation Basis Algorithm 51 ) Used in chunk 50. 
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